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Discrete complex analysis on planar quad-graphs 


Alexander I. Bobenko*’^ 
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Abstract 

We develop a linear theory of discrete complex analysis on general quad-graphs, continuing and 
extending previous work of DufSn, Mercat, Kenyon, Chelkak and Smirnov on discrete complex 
analysis on rhombic quad-graphs. Our approach based on the medial graph yields more instructive 
proofs of discrete analogs of several classical theorems and even new results. We provide discrete 
counterparts of fundamental concepts in complex analysis such as holomorphic functions, deriva¬ 
tives, the Laplacian, and exterior calculus. Also, we discuss discrete versions of important basic 
theorems such as Green’s identities and Cauchy’s integral formulae. For the first time, we discretize 
Green’s first identity and Cauchy’s integral formula for the derivative of a holomorphic function. 
In this paper, we focus on planar quad-graphs, but we would like to mention that many notions 
and theorems can be adapted to discrete Riemann surfaces in a straightforward way. 

In the case of planar parallelogram-graphs with bounded interior angles and bounded ratio of side 
lengths, we construct a discrete Green’s function and discrete Cauchy’s kernels with asymptotics 
comparable to the smooth case. Further restricting to the integer lattice of a two-dimensional 
skew coordinate system yields appropriate discrete Cauchy’s integral formulae for higher order 
derivatives. 
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1 Introduction 

We start with a short history of discrete complex analysis before sketching the organization of our 
paper. Our aim is to give just a very rough overview of the linear discrete theories of complex analysis 
we are enhancing and to name connections to statistical physics and the nonlinear discrete theory of 
complex analysis based on circle patterns. For a more detailed discussion of these topics, we refer to 
the survey of Smirnov [23] , on which our introduction is based and that also discusses applications to 
probability theory and mathematical physics, and to the book [2| of the first author and Suris that not 
only investigates the influence of discrete integrable systems to discrete theories of complex analysis, 
but also gives a good overview of discrete differential geometry in general. 

Linear theories of discrete complex analysis look back on a long and varied history. Already Kirch- 
hoff’s circuit laws describe a discrete harmonicity condition for the potential function whose gradient 
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describes the current flowing through the electric network. A notable application of Kirchhoff’s laws in 
geometry was the article [3] of Brooks, Smith, Stone, and Tutte, who used coupled discrete harmonic 
functions (in fact, discrete holomorphic functions) to construct tilings of rectangles into squares with 
different integral side lengths. 

Discrete harmonic functions on the square lattice were studied by a number of authors in the 1920s, in¬ 
cluding Courant, Friedrichs, and Lewy, who showed convergence of solutions of the Dirichlet boundary 
value problem to their corresponding continuous counterpart [?]• Recently, Skopenkov studied Dirich¬ 
let boundary value problems of discrete analytic functions on general quad-graphs in [22], and he 
showed a convergence result in the case of quad-graphs where the diagonals of quadrilaterals intersect 
orthogonally. 

Discrete holomorphic functions on the square lattice were studied by Isaacs m- He proposed two 
different definitions for holomorphicity. The first one said that / is monodiffric of the first kind if 

f{z + is) - f{z) = i{f{z + e)- f{z)), 

where s denotes the side length of the squares. This definition is not symmetric on the square lattice, 
but it becomes symmetric on the triangular lattice obtained by inserting the diagonals parallel to 1 — i. 
Dynnikov and Novikov studied this notion in |10] . 

Isaac’s second definition was given by 

f{z + is) - f{z -I- e) = i{f{z -h (1 -h i)e) - f{z)). 

In this context, it is natural to consider the real part of a discrete holomorphic function as being 
defined on one type of vertices, say black, and the imaginary part on the other type of vertices, say 
white, corresponding to a bipartite decomposition of the square lattice. Lelong-Ferrand reintroduced 
this notion in m and developed the theory to a level that allowed her to prove the Riemann mapping 
theorem using discrete methods m- Duffin also studied discrete complex analysis on the square grid 
[8] and was the first who extended the theory to rhombic lattices |9|. Mercat |18] . Kenyon |15] . and 
Chelkak and Smirnov [5] resumed the investigation of discrete complex analysis on rhombic lattices 
or, equivalently, isoradial graphs. In these settings, it was natural to split the real and the imaginary 
part of a discrete holomorphic function to the two vertex sets of a bipartite decomposition. 

Some two-dimensional discrete models in statistical physics exhibit conformally invariant properties in 
the thermodynamical limit. Such conformally invariant properties were established by Smirnov for site 
percolation on a triangular grid [23| and for the random cluster model |25|, by Chelkak and Smirnov 
for the Ising model |6], and by Kenyon for the dimer model on a square grid (domino tiling) [T3|. In 
all cases, linear theories of discrete analytic functions on regular grids were highly important. Kenyon 
as well as Chelkak and Smirnov obtained important analytic results mu that were instrumental in 
the proof that the critical Ising model is universal, i.e., that the scaling limit is independent of the 
shape of the lattice [6]. Already Mercat related the theory of discrete complex analysis to the Ising 
model and investigated criticality [TH] . 

Important non-linear discrete theories of complex analysis involve circle packings or, more generally, 
circle patterns. Rodin and Sullivan first proved that the Riemann mapping of a complex domain 
to the unit disk can be approximated by circle packings m- A similar result for isoradial circle 
patterns, even with irregular combinatorics, is due to Bucking [Ij. In their paper [T], the hrst author, 
Mercat, and Suris showed how the linear theory of discrete holomorphic functions on quad-graphs 
can be obtained by linearizing the theory on circle patterns: Discrete holomorphic functions describe 
infinitesimal deformations of circle patterns. In the case of parallelogram-graphs, they embedded the 
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quad-graph in Z” and introduced a discrete exponential function and a discrete logarithm, generalizing 
Kenyon’s discrete exponential and discrete Green’s function HU. 

Our setup in Section [2] is a strongly regular cellular decomposition of the complex plane into quadri¬ 
laterals, called quad-graph, that we assume to be bipartite. Of crucial importance for our work is the 
medial graph of a quad-graph. It provides the connection between the notions of discrete derivatives 
of Kenyon m, Mercat [19], and Chelkak and Smirnov [5], extended from rhombic to general quad- 
graphs, and discrete differential forms and discrete exterior calculus as suggested by Mercat [THl 1^ . 
Concerning discrete differential forms, we get essentially the same definitions as Mercat proposed in 
|20j . However, our notation of discrete exterior calculus is slightly more general and shows its power 
when considering integral formulae. Discrete Stokes’ Theorem 12.8.11 will be a consequence and not 
part of the definition of the discrete exterior derivate as in the work of Mercat [181 [20] , and in Theo¬ 
rem [ 23 S] we prove that the discrete exterior derivative is a derivation of the discrete wedge-product. 
These two theorems are the most powerful ones in our setting. In particular, the proof of discrete 
Green’s identities in Theorem 12.4.41 is an immediate corollary. Here, a discrete version of Green’s first 
identity is provided for the first time. 

According to the paper |5] of Chelkak and Smirnov, one of the unpleasant facts of all discrete theories 
of complex analysis is that (pointwise) multiplication of discrete holomorphic functions does not yield 
a discrete holomorphic function in general. We suggest at least a partial solution in Corollary 12.8.51 
where we describe how the medial graph allows to (kind of pointwise) multiply discrete holomorphic 
functions to a function that is not defined on the vertex sets of the original graphs anymore, but that 
is discrete holomorphic in a certain sense. 

Based on Skopenkov’s results on the existence and uniqueness of solutions to the discrete Dirichlet 
boundary value problem [22], we prove surjectivity of the discrete differentials and the discrete Lapla- 
cian seen as linear operators in Theorem 12. 4. Ill In particular, discrete Green’s functions and discrete 
Cauchy’s kernels z~^ exist. As a consequence, we can formulate discrete Cauchy’s integral formu¬ 
lae for discrete holomorphic functions in Theorem 12.6.21 and for the discrete derivative of a discrete 
holomorphic function on the vertices of the quad-graph in Theorem 12.6.81 Unfortunately, we cannot 
require any certain asymptotic behavior of them in the general setup so far. But at least we show in 
Theorem 12.4.121 that any discrete harmonic function with asymptotics is constant, provided 

that all interior angles and side lengths of the quadrilaterals are bounded. 

Section |3| is devoted to discrete complex analysis on planar parallelogram-graphs. There, we construct 
discrete Green’s functions and discrete Cauchy’s kernels with asymptotics similar to the functions in 
the rhombic case |15[ H] [5] and close to the smooth case. This will be proven in Theorems 18.2.21 
18.8.11 and 18.8.81 assuming that the interior angles and the ratio of side lengths of all parallelograms 
are bounded. The construction of these functions is based on the discrete exponential introduced by 
Kenyon on quasicrystallic rhombic quad-graphs |15j and by the first author, Mercat, and Suris on 
quasicrystallic parallelogram-graphs [T]. In the end, we close with the very special case of the integer 
lattice of a skew coordinate system in the complex plane. In this case, we show in Theorem 18.4.11 
that discrete Cauchy’s integral formulae for higher order discrete derivatives of a discrete holomorphic 
function exist and that the asymptotics of these formulae match the expectations from the previous 
results. 

Finally, in the appendix we give the proofs of some statements concerning the combinatorics of 
parallelogram-graphs that we use for determining the asymptotic behavior of certain discrete functions. 
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2 Discrete complex analysis on planar quad-graphs 


Although we focus on planar quad-graphs in this paper, many of our notions and theorems generalize 
to discrete Riemann surfaces. A corresponding linear theory of discrete Riemann surfaces will be 
discussed in a subsequent paper and can also be found in the PhD thesis of the second author [12]. 

2.1 Basic definitions and notation 

The aim of this section is to introduce first bipartite quad-graphs and some basic notation in Sec¬ 
tion 12.1.11 and then to discuss the medial graph in Section 12.1.21 


2.1.1 Bipartite quad-graphs 


We consider a strongly regular cellular decomposition of the complex plane C described by an embed¬ 
ded bipartite quad-graph A such that 0-cells correspond to vertices R(A), 1-cells to edges E{A), and 
2-cells to quadrilateral faces F{A). We refer to the maximal independent sets of vertices of A as black 
and white vertices. Furthermore, we restrict to locally finite cellular decompositions, i.e., a compact 
subset of C contains only finitely many quadrilaterals. 

The assumption of strong regularity asserts that the boundary of a quadrilateral contains a particular 
vertex or a particular edge at most once and that two different edges or faces have at most one vertex 
or edge in common, respectively. As a consequence, if a line segment connecting two vertices of A is 
the (possibly outer) diagonal of a quadrilateral, there is just one such face of A. Let T and T* be the 
graphs defined on the black and white vertices where the edges are exactly the corresponding diagonals 
of faces of A. If the diagonal lies outside the face, it is more convenient to consider the corresponding 
edge of T or T* not as the straight line segment connecting its two endpoints, but as a curve lying 
inside the face. Then, the duality between T and T* becomes obvious: Black and white vertices are 
in one-to-one correspondence to the white and black faces they are contained in, and black and white 
edges dual to each other are exactly these who cross each other. For simplicity, we identify vertices of 
A or of F and F* with their corresponding complex values, and to oriented edges of A, F, F* we assign 
the complex numbers determined by the difference of their two endpoints. 

To A we associate its dual 0 := A*. In general, we look at 0 in an abstract way, identifying vertices 
or faces of 0 with corresponding faces or vertices of A, respectively. However, in the particular case 
that all quadrilaterals are parallelograms, it makes sense to place the vertices of <0 at the centers of 
the parallelograms. Here, the center of a parallelogram is the point of intersection of its two diagonals. 
Further details will be given in Sections 12.2.II and 12.2.31 

If a vertex v E R(A) is a vertex of a quadrilateral Q E F{A) = R(0), we write Q ^ v oi v ^ Q 
and say that v and Q are incident to each other. The vertices of Q are denoted by , re-, 6+, t(;+ in 
counterclockwise order, where h± E R(r) and w± E F(r*). 

Definition. For a quadrilateral Q E R(0) we define 


Moreover, let 



ipQ := arccos 




= arccos 



( ib+ - b-){w+ - W-) \\ 
\ \b+ -b-\\w+ -W-\ j j 
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be the angle under which the diagonal lines of Q intersect. 

Note that 0 < (pg < vr. Figure [T] shows a finite bipartite quad-graph together with the notations we 
have introduced for a single quadrilateral Q and the notations we are using later for the star of a 
vertex v, i.e., the set of all faces incident to v. 


V2 



Figure 1: Bipartite quad-graph with notations 


In addition, we denote by Oo always a connected subset of 0- It is said to be simply-connected if the 
corresponding set of cells in C is simply-connected. Its vertices induce a subgraph Aq of A together 
with subgraphs Fq of F and Fq of F*. For simplicity, we assume that the induced subgraphs are 
connected as well. By SAq we denote the subgraph of Aq that consists of boundary vertices and edges. 

2.1.2 Medial graph 

Definition. The medial graph X of A is defined as follows. Its vertex set is given by all the midpoints 
of edges of A, and two vertices are adjacent if and only if the corresponding edges belong to the same 
face and have a vertex in common. The set of faces of X is in bijection with F(A) U 1^(0): A face 
Fy corresponding to n G I^(A) consists of the midpoints of edges of A incident to n, and a face Fq 
corresponding to Q G 1^(0) consists of the midpoints of the four edges of A belonging to Q. 

Any edge e of A" is the common edge of two faces Fq and for Q ~ n, denoted by [Q, n]. 

Let Q (0) and vq ~ Q. Due to Varignon’s theorem, Fq is a parallelogram and the complex number 
assigned to the edge e = [Q,no] connecting the midpoints of edges vqv'_ and vqv'j^ of A is just half of 
e = v'_^_ — v'_. Note that if Q is nonconvex, some part of Fq lies outside Q and it may happen that v 
lies not inside Fy. These situations do not occur if all quadrilaterals are convex. 

In Figure [21 showing A with its medial graph, the vertices of Fq and Fy are colored gray. 

For a subgraph Oo F 0, we denote by Xq C X the subgraph of X consisting of all edges [(5,n] where 
Q G IL(Oo) and v ^ Q. In the case of convex quadrilaterals, this means that we take all edges of X 
lying inside Oo- The medial graph X corresponds to a (strongly regular and locally finite) cellular 
decomposition of C in a canonical way. In particular, we can talk about a topological disk in F{X) 
and about a (counterclockwise oriented) boundary dX^. 

Definition. For v G IL(A) and Q G D(0), let Py and Pq be the closed paths on X connecting the 
midpoints of edges of A incident to v and Q, respectively, in counterclockwise direction. In Figure El 
their vertices are colored gray. We say that Py and Pq are discrete elementary cycles. 
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Figure 2: Bipartite quad-graph (dashed) with medial graph (solid) 


2.2 Discrete holomorphicity 

In the classical theory, a real differentiable function f : U C C —C is holomorphic if and only if the 
Cauchy-Riemann equation dxf = —idyf is satisfied in all points of U. Moreover, holomorphic functions 
with nowhere-vanishing derivative preserve angles, and at a single point, lengths are uniformly scaled. 



Figure 3: Discretization of Cauchy-Riemann equation 


Now, let us imagine a rhombus Q in C with vertices b-,w-,b+,w+ and diagonals b-b^,w-w^ aligned 
to the real x- and the imaginary y-axis, see Figure [3l Then, 

discretizes dxf{Q), 

b+ — b- 

fiw+) - fiw-) _iQ 

W-^- — W- 

This motivates the following definition of discrete holomorphicity due to Mercat [20] that was also 
used previously in the rhombic setting by Duffin [S] and others. 

Definition. Let Q G R(0) and / be a complex function on b-,w-,b+,'w+. f is said to be discrete 
holomorphic at Q if the discrete Cauchy-Riemann equation is satisfied: 

f{b+) - f{b-) ^ f{w+) - f{w-) 

6+ — 6_ w+ — w- 

If Q is not a rhombus, we do not have an immediate interpretation of the discrete Cauchy-Riemann 
equation as before. But if a discrete holomorphic function / does not have the same value on both 
black vertices and 6+, it preserves the angle ipg and / uniformly scales the lengths of the diagonals 
of Q. However, the image of Q under / might be a degenerate quadrilateral. 
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We immediately see that for discrete holomorphicity, only the differences at black and at white vertices 
matter. Hence, we should not consider constants on H(A), but biconstants [19] determined by each a 
value on H(r) and H(r*). A function that is constant on H(r) and constant on H(r*) is said to be 
biconstant. 

In the following, we will define discrete analogs of d, d, first of functions on V (A) in Section 12.2.11 
and later of functions on V{()) in Section 12.2.31 Before, we introduce discrete differential forms in 
Section 12.2.21 


2.2.1 Discrete derivatives of functions on H(A) 


Remember that the derivatives d := dz and d := dz in complex analysis are defined through 


d = {dx — idy) /2 and d = {dx + idy) /2; 


the coefficients in front of the partial derivatives in d are complex conjugate to the coefficients ap¬ 
pearing in d. Using the interpretation 


f{b+)-f{b.) 

b+-b- 

f{w+) - fjw-) 

— W- 


dxf, 


-idyf 


for rhombi Q as above, the definition of discrete derivatives in rhombic quad-graphs used by Chelkak 
and Smirnov in |5|, 


OaHQ) 

4/(Q) 


1 

2 

1 

2 


f{b+)-f{b-) ^ f{w+) - f{w-) \ 

6 + — 5 _ — w- J 

/( 6 +)-/( 6 -) ^ f{w+)- f{w-) \ 

b+ - b- w+ - w- J 


becomes plausible. 

In a general quadrilateral Q, the diagonals usually do not intersect orthogonally and their quotient 
is not necessarily purely imaginary. Therefore, we have to choose different factors in front of the 
difference quotients, taking the deviation ipq — 7r/2 from orthogonality into account. 

Definition. Let Q € U(0), and let / be a complex function on 5_, 6+, u;+. The discrete derivatives 

9a/, 9a/ are defined by 


9a/(Q) := 
BaHQ) 


X f{b+)-f{b-) , ^ f{w+)-f{w-) 

f{w+) - fiw-) 


b+-b. 

, /(M-/(fe-) , . 

- , , -+ 

0-1- — b- 


— W- 


where 2 Aq := exp {-i [ifg - |)) / sin((/9Q). 


Obviously, biconstant functions have vanishing discrete derivatives. If the quadrilateral Q is a rhombus, 
ipQ = 7r/2 and Xq = 1/2. Thus, the definition above then reduces to the previous one in |5]. The 
definition of discrete derivatives also matches the notion of discrete holomorphicity; and the discrete 
derivatives approximate their smooth counterparts correctly up to order one for general quad-graphs 
and up to order two for parallelogram-graphs: 
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Proposition 2.2.1. Let Q G P(0) and f be a complex function on b-,w-,b^,w+. 

(i) f is discrete holomorphic at Q if and only if d\f{Q) = 0. 

(a) If f{v) = V, dAf{Q) = 0 and 5a/(Q) = 1. 

(in) If Q is a parallelogram and f{v) = , dxfiQ) = 0 and d^fiQ) = “IQ, where Q as a vertex is 

placed at the center of the parallelogram. 

(iv) If Q is a parallelogram and f{v) = \v\^, 5a/(Q) = d\f{Q) = Q, where Q as a vertex is placed 
at the center of the parallelogram. 

Proof, (i) We observe that 

2sin((^Q)AQ ^ exp(i (v?q - f)) ^ w+ - W- 

b+-b- b+-b- \w+- w-\\b+- b-\' 

2sm{ipQ)XQ ^ exp(-i {ipg - f)) ^ b+ - b- 

w+— W- w+— W- \w+— W-\\b+— b-\ 

So if we multiply d^fiQ) by 2i\w+ — rr_||6+ — 6_| sin(99Q) / 0, we obtain 

(u;+ - W-) (/ (6+) - / (6_)) - (6+ - 6_) (/ (u;+) - / («;_)). 

The last expression vanishes if and only if the discrete Cauchy-Riemann equation is satisfied. 

(ii) Clearly, f{v) = v satisfies the discrete Cauchy-Riemann equation. By the first part, d\f{Q) = 0. 
Due to 2sin((/9Q) = exp {-i {(pq - f)) + exp [i [pq - f)), 9 a/(( 5) simplifies to Aq -h Aq = 1 . 

(hi) For f{v) = the discrete Cauchy-Riemann equation is equivalent to 6+ -|- 6_ = + W-. But 

since Q is a parallelogram, both {b+ + b-)/2 and {w+ + W-)/2 equal the center Q of the parallelogram. 
Thus, / is discrete holomorphic at Q and 

dAf (Q) = I^Qib+ + b-) + Xq{w+ + W-) = 2Q{Xq + Xq) = 2Q. 

(iv) Since / is a real function, dAfiQ) = <9a/(Q) follows straight from the definition. Let z G C be 
arbitrary. If g{v) := vz, dAg{Q) = ^ and dAg{Q) = 0 by the second part. So if we define the function 
h{v) := |u — zp = |up — vz — vz + \z\‘^, then dAh{Q) = cIa/(Q) — Hence, the statement is invariant 
under translation, and it suffices to consider the case when the center of the parallelogram Q is the 
origin. Then, 6+ = —6_ and w+ = —W- since Q is a parallelogram. It follows that f{b-) = fipjf) 
and f{w-) = f{w+), so dAf{Q) = 0 as desired. □ 

Our first discrete analogs of classical theorems are immediate consequences of the discrete Cauchy- 
Riemann equation; 

Proposition 2.2.2. Let Oo ^ 0 and f : V{A) —)■ C 6e discrete holomorphic. 


(i) If f is purely imaginary or purely real, f is biconstant. 

(ii) If dAf = 0, / is biconstant. 



Proof, (i) Let 6 _, 6 + be two adjacent vertices of L, and let b-,'W-,b+,w+ E l^(Ao) be the vertices 
of the corresponding quadrilateral in 0- Let us assume that /( 6 +) / fib-). Due to the discrete 
Cauchy-Riemann equation, 

fiw+) - fjw-) ^ w+-w_ 

/( 6 +) - fib.) 6 + - 6 _ • 

The left hand side is real, but the right hand side is not, contradiction. The same argumentation goes 
through if we start with two adjacent vertices of T*. 

(ii) Since / is discrete holomorphic, 

fib+)-fib.) fiw+)-fiw.) _ 
b+ — b- w-\- — w- 

d\f = 0 then yields that both sides of the discrete Cauchy-Riemann equation equal zero, so / is 
constant on R(ro) and constant on R(rQ). □ 

2.2.2 Discrete differential forms 

We mainly consider two type of functions, functions / : R(A) —C and functions h : R(<C>) —> C. An 
example for a function on the quadrilateral faces is d{^f. 

A discrete one-form w is a complex function on the oriented edges of the medial graph X, and a 
discrete two-form D is a complex function on the faces of X. The evaluations of uj at an oriented edge 
e of W and of D at a face F of X are denoted by oj and D, respectively. 

If P is a directed path in X consisting of oriented edges ei, 62 ,..., e„, the discrete integral along P is 
defined as fpCO = le^ closed paths P, we write fpCO instead. In the case that P is the 

boundary of an oriented disk in X, we say that the discrete integral is a diserete contour integral with 
discrete contour P. Similarly, the discrete integral of D over a set of faces of X is defined. 

It turns out that discrete one-forms that actually come from discrete one-forms on T and T* are of 
particular interest. We say that such a discrete one-form oj is of type (). It is characterized by the 
property that for any Q E R(0) there exist complex numbers p,q such that lo = pdz + qdz on all 
edges e = [Q,v], V ^ Q. In an analogous fashion we can define discrete one-forms u' of type A by 
representing them as u' = pdz + qdz on all edges e = [n, Q], Q ~ n, for some fixed vertex u E R(A) and 
complex numbers p, q attached to it. However, these discrete one-forms do not play such an important 
role than the others. We mainly give them to get a complete picture. 

As we are mainly interested in functions / : R(A) —>■ C and h : R(0) —^ C, discrete two-forms of 
particular interest are those that vanish on faces of X corresponding to vertices of either 0 or A. 
We will call them of being of type A or type (), respectively. These discrete two-forms correspond to 
functions on R(A) or R(0) by the discrete Hodge star that will be defined later in Section [2.,‘I.,'II 

Definition. Let P be a face of the medial graph X. We define ar(P) to be twice the Euclidean area 
of P in the complex plane. In contrast, area(P) will always denote the Euclidean area of a polygon 

P. 

Remark. As we have mentioned before, our main objects either live on the quad-graph A or on its 
dual 0- Thus, we have to deal with two different cellular decompositions at the same time. The medial 
graph has the crucial property that its faces are in one-to-one correspondence to vertices of A and of 
0, i.e., to faces of 0 and of A. Furthermore, the Euclidean area of the Varignon parallelogram inside 
Q G V (0) is just half of the area of Q. In some sense, a corresponding statement is true for the cells 
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of X corresponding to vertices of A, i.e., faces of 0. However, there is no natural embedding of 0 in 
the general setting. But in the particular case of parallelogram-graphs, we can make the statement 
precise: If an edge QQ' of 0 is represented by the two line segments that connect the centers of the 
parallelogram Q and Q' with the midpoint of their common edge, then the Euclidean area of the face 
of X corresponding to a vertex v £V (A) = F(0) is exactly half of the area of the face v. 

In summary, the medial graph allows us to deal with just one decomposition of the complex plane, 
but we have to count areas twice in order to get the right prefactors as in the continuous setup. 


For later purposes, let us compute ar(E’). First, if F is the Varignon parallelogram corresponding to 
the quadrilateral Q E F(0), 

ar(^) = ^\b+- b_\\w+ - u;_| sin(v;Q). 

Note again that the Euclidean area of F is just half of ar (F). Second, if F is the face of X corresponding 
to a vertex v G V (A), its Euclidean area is a quarter of the area of the polygon v'iv '2 ... in the star 
of V. So ar(F) equals 

i ^ area(A'(;n'_iu') = ^ ^ Im (^(u' - u) ^ X] > 

Qa'^v Qs''^v Qs~'y 

using that + vv'^) = Eg,-,, (vF, + vv',) is real. 

Definition. The discrete one-forms dz and dz are given by f^dz = e and J^dz = e for any oriented 
edge e of X. The discrete two-forms Ha and are zero on faces of X corresponding to vertices of 0 
or A, respectively, and defined by 


Ha = —2iar(F) and 


Ha = 


JJF JJF 

on faces F corresponding to vertices of A or 0, respectively. 
Euclidean area of F. 


—2iar(F) 

As defined above, ar(F) is twice the 


Remark. Ha and H^ are the straightforward discretizations of dz A dz having in mind that they are 
essentially defined on faces of 0 or of A, respectively. It turns out that in local coordinates, we can 
perform our calculations with Ha and in the discrete setting exactly as we do with dz A dz in the 
smooth theory. Indeed, we will see in Section 12.3.21 that Ho is indeed the discrete wedge product of 
dz and dz seen as discrete one-forms of type 0- The same is true for Ha if we consider dz and dz as 
being of type A, but the discrete wedge-product is of interest just for discrete one-forms of type 0- 

Definition. Let / : E(A) —C, /i : E(0) ^ C, w a discrete one-form, and Hi,H 2 discrete two-forms 
of type A and 0, respectively. For any edge e = [Q,v] and any faces Fy,FQ of X corresponding to the 
vertex star of u € E(A) or the Varignon parallelogram inside Q E E(0), we define the products fuj, 
huj, /Hi, and /iH 2 by 


JF: 

= f{v) j UJ 

and 

II 

■■= f{v] 

'II^^- 

II 

:=0; 

e 

e 


Fy 


Fy 

Fq 


h- 

= h{Q) j uj 

and 


:=0, 


IIhf- 

= h{Q) 11^2 

e 

e 


Fy 



Fq 

Fq 
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In the following table, we give a quick overview of various discrete differential forms most of them will 
be discussed in Section 12.31 and state whether they are essentially functions on V (A) or functions on 
y(0) or entirely objects on the cellular decomposition X. Although discrete one-forms of type A or of 
type 0 do not live themselves on A or 0, they are described by two functions defined on the vertices 
of A or 0, respectively. 



A 

0 

X 

0-forms 

f,g:V{A)^€ 

hi,h2:V{0)^€ 

f -9 = Jifdg + gdf) 


d()h, dg,h 

dkf, dxf 


1-forms 

dh 

df 

fdg + gdf 


hidz + h2dz 

fdz+gdz 

fhdz 


r] of type A 

a;, a;' of type 0 

foJ 

2-forms 


Qo 



*/ 

'kh 



du 

dr] 

d{fhdz) 


fdu 

CO Alo' 

d{fuj) 


2.2.3 Discrete derivatives of functions on 1^(0) 

Before we pass on to discrete derivatives of functions on V(0), we first prove an alternative formula 
for the discrete derivatives of functions on V (A). 

Lemma 2.2.3. Let Q E 1^(0) and f be a complex function on the vertices of Q. Let Pq be the discrete 
elementary cycle around Q and F the face of X corresponding to Q. Then, 

Pq 

Pq 

Proof. Since F is a parallelogram, /(6+) and —/(6_) are multiplied by the same factor (tt)+ — W-)f2 
when evaluating the discrete contour integral fdz. Therefore, the prefactor in front of /(5+)—/(6_) 
in the right hand side of the first equation is 

. W+ - W- _ _ w+ - W- __ exp [-i [ipQ - D) _ Aq 

4ar(F) 2sin((/?Q)|t(;+- t(;_||6+- 6_| 2sin((/7Q)(6+- 6_) 6+- 6_ 
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(compare with the proof of Proposition l2.2.I] (i)). which is exactly the prefactor appearing in d\f{Q). 
In an analogous manner, the prefactors in front of /(rc+) — fiw-) are equal. This shows the first 
equation. The second one follows from the first, noting that the prefactors in front of /(6+) —/(6_) and 
f{w+) — f{w-) on both sides of the second equation are just complex conjugates of the corresponding 
prefactors appearing in the first equation. □ 


w+ 




(a) Lemma 12.2.31 for 9a, Sa (b) Definitions of 9o,9o 

Figure 4; Integration formulae for discrete derivatives 


Inspired by Lemma 12.2.31 we can now define the discrete derivatives for complex functions on 14(0), 
see Figure 0] (b). 

Definition. Let v G 14(A) and /i be a complex function defined on all quadrilaterals Qs ~ v. Let 
be the discrete elementary cycle around v and F the face of X corresponding to v. Then, the discrete 
derivatives doh, d(^h aX v are defined by 

Pv 

Pv 


h is said to be discrete holomorphic at v if d(^h{v) = 0. 

Note that in the rhombic case, our definition coincides with the one used by Chelkak and Smirnov in 
[5]. As an immediate consequence of the definition, we obtain a discrete Morera’s theorem. 

Proposition 2.2.4. Functions f : 14(A) — )• C and h : 14(0) — >■ C are discrete holomorphic if and only 
if §p fdz = 0 and §p hdz = 0 for all discrete contours P. 

Proof. Clearly, fp fdz = f{v) ^p dz = 0 for any discrete elementary cycle around a vertex v of 
14(A). Similarly, fp^ hdz = 0 for any Q G 14(0). Using Lemma r2 .2.3 1 and the definition of 5o, / and h 
are discrete holomorphic if and only if fp fdz = 0 and fp hdz = 0 for all discrete elementary cycles P. 
To conclude the proof, we observe that any integration along a discrete contour can be decomposed 
into integrations along discrete elementary cycles. □ 

The discrete derivatives of constant functions on 14(0) vanish. As an analog of Proposition 12.2.11 we 
prove that the discrete derivatives <9o,9<) locally approximate their smooth counterparts correctly up 
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to order one if the vertices Q are placed at the midpoints of black or white edges. Note that even for 
rhombic quad-graphs, these discrete derivatives generally do not coincide with the smooth derivatives 
in order two. 

Proposition 2.2.5. Let v € 1^(A), and let h he a complex function on all faces Qg ~ v. Assume that 
the vertices Qg are placed at the midpoints (Ws_i +Vg)/2 of edges of the star of v. Then, if h{Q) = Q, 
d()h{v) = 0 and d(^h{v) = 1. 

Proof. 


4:j)hdz= - v'g-i) = (K)^ - K-i)^) = 0, 

p Qa'^V Qs'^V 



= —2i I™ = —8iai{Fy). 

Qa'-^V 


These equations yield d()h{v) = 0 and d()h{v) = 1. □ 

Remark. As before, parallelogram-graphs play a special role. In a parallelogram-graph, the midpoint 
of a black edge equals the midpoint of its dual white edge. Placing the vertices Q £ V (0) at the centers 
of the corresponding parallelograms then yields a global approximation statement. 


In [5], Chelkak and Smirnov used averaging operators to map functions on I^(A) to functions on P(0) 
and vice versa. On parallelogram-graphs, the averaging operator m{f){Q) := Ylvr^Q f ^ actually 
maps discrete holomorphic functions / : I^(A) —>■ C to discrete holomorphic functions on 1^(0)- The 
corresponding statement for rhombic quad-graphs was shown in [S], our proof is similar. 

Proposition 2.2.6. Let A be a parallelogram-graph and f : P(A) -£■ C be discrete holomorphic. Then, 
m{f) : P(0) ^ C is discrete holomorphic. 


Proof. Let us consider the star of the vertex v £ P(A). Since / is discrete holomorphic, the discrete 
Cauchy-Riemann equation is satisfied on any Qg ~ v. Therefore, we can express f{vg) in terms of 
f{v), f{Vg) and f{Vg_i). Plugging this in the definition of the averaging operator, we obtain 


4m(/)(Q,) = 2f{v) +-+- ffLZTffr - 


V' — v 


s-1 


s—1 


V' — V 


Vg-l - V 


= 2f{v) + 2 ; f{v',) - 2 ^ /K-i). 


V' — v 


s—1 


V' — V 


s-1 


Here, we have used the properties Vg — v'g_i = Vg — v and Vg — v'g = v'g_i — u of the parallelogram Qg. 
Therefore, 


2 


m{f)dz 


Pv 



-v)fiv'g) 


Y (^^-1 “ ^)/ K - i ) = 

Qa'^H 


and m{f) is discrete holomorphic at v. 


□ 
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Remark. As noted by Chelkak and Smirnov, discrete holomorphic functions on V{()) cannot be 
averaged to discrete holomorphic functions on V (A) in general. 


As mentioned above, our main interest lies in functions that are defined either on the vertices or the 
faces of the quad-graph. Now, extending / ; 1^(A) —>■ C to a complex function on F{X) by using 
its average m{f) on V{()) seems to be an option. However, functions on H(A) and on H(0) behave 
differently. In Corollary 12.3.31 we will see that d\f is discrete holomorphic if / is, but does 

not need to be discrete holomorphic in general. So to make sense of differentiating twice, we can only 
consider functions on H(A). 

Definition. Let /i ,/2 : H(A) —C and : H(0) ^ C. Their discrete scalar products are defined 

as 

(/i,/ 2 ) := J j fih^K and (/ii,/i 2 ) := jj hih 2 ^(^, 

F(X) Fix) 

whenever the right hand side converges absolutely. 

Proposition 2.2.7. —d<) and —d<) are the formal adjoints of Ba and d\, respectively. That is, if 
f : H(A) —)• C or /i : H(0) —C is compactly supported, 

{d\f, h) + (/, Bf^h) = 0 = (Sa/, h) + (/, Boh). 

Proof. Using Lemma [2.2.31 and B^h = B^h, we get 

—2i{BAf,h) — 2i{f,Boh)= ^ h{Q)<lfdz+ ^ f{v)([hdz=([fhdz = 0, 

QeViO) ^ veViA) i( i 

where P is a large contour enclosing all the vertices of A and 0 where / or /i do not vanish. In 
particular, fh vanishes in a neighborhood of P. In the same way, {B^f, h) -|- (/, B^h) = 0. □ 

Remark. In their work on discrete complex analysis on rhombic quad-graphs, Kenyon m and Mercat 
m defined the discrete derivatives for functions on the faces in such a way that they were the formal 
adjoints of the discrete derivatives for functions on the vertices of the quad-graph. 


In Corollary 12.3.31 we will prove that B^f is discrete holomorphic if / ; U(A) —)• C is. Conversely, 
we can find discrete primitives of discrete holomorphic functions on simply-connected domains Ooj 
extending the corresponding result for rhombic quad-graphs given in the paper of Chelkak and Smirnov 


0 . 


Proposition 2.2.8. Let Oo F () be simply-connected. Then, for any discrete holomorphic function h 
on U(<0o); there is a discrete primitive f := J h on U(Ao), i.e., f is discrete holomorphic and B\f = h. 
f is unique up to two additive constants on Lq and Lq. 


Proof. Since h is discrete holomorphic, <fp hdz = 0 for any discrete contour P. Thus, hdz can be 
integrated to a well defined function fx on V{X) that is unique up to an additive constant. Using 
that hdz is a discrete one-form of type <0, we can construct a function f on V (A) such that the equation 
fx {{v + w) /2) = (/(u) -t- f{w)) /2 holds for any edge {v,w) of A. Given fx, f is unique up to an 
additive constant. 
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In summary, / is unique up to two additive constants that can be chosen independently on Fq and 
Fq. By construction, / satisfies 

f{b+)-f{b-) ^ ^ f{w+)-f{w_) 

6+ — b- w-\- — ic¬ 
on any quadrilateral Q. So / is discrete holomorphic and d\f = h. □ 

2.3 Discrete exterior calculus 

Our notation of discrete exterior calculus is similar to the approach of Mercat in [T8l [T9l [20], but 
differs in some aspects. The main differences are due to our different notation of multiplication of 
functions with discrete one-forms, which allows us to define a discrete exterior derivative on a larger 
class of discrete one-forms in Section 12.3.11 It coincides with Mercat’s discrete exterior derivative in 
the case of discrete one-forms of type 0 that Mercat considers. In contrast, our definitions are based 
on a coordinate representation, making the connection to the smooth case evident. Eventually, they 
lead to essentially the same definitions of a discrete wedge product in Section 12.3.21 and a discrete 
Hodge star in Section 12.3.31 as in [20] . 

2.3.1 Discrete exterior derivative 

Definition. Let / : E(A) ^ C and h : E(<}) —^ C. We define the discrete exterior derivatives df and 
dh as follows: 


df := dxfdz + dAfdz; 
dh := d()hdz + d(^hdz. 

Let a; be a discrete one-form. Around faces and Fq of X corresponding to vertices v G V (A) and 
Q G V (0), respectively, we write tv = pdz + qdz with functions p, q defined on faces ~ u or vertices 
b±,w± Q, respectively. The discrete exterior derivative doj is given by 

duj\F^ := {d^q - d^p) Ha, 
duj\FQ := {dAq - Bap) 

The representation of iv as pdz+qdz (p, q defined on edges of X) we have used above may be nonunique. 
However, duj is well defined by discrete Stokes’ theorem that also justifies our definition of df and div. 
Note that Mercat defined the discrete exterior derivative by the discrete Stokes’ theorem [18] . 

Theorem 2.3.1. Let f : E(A) — )• C and tv be a discrete one-form. Then, for any directed edge e of 
X starting in the midpoint of the edge vv'_ and ending in the midpoint of the edge vv'_^_ of A and for 
any face F of X with counterclockwise oriented boundary dF we have: 



fK) - fjv'-) 

2 


f{v) + fiv+) _ f{v) + f{v-) _ 

2 2 



F 


du 



dF 


VJ. 
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Proof. Let V- be the other vertex of the quadrilateral Q with vertices v, v'_ and v'j^. Without loss of 
generality, let v be white. Then, df equals 




v'-v' 


+ 9a/ 


v\ — v' 


+ ■^q)(/(^+) “ fW-)) + 2 ( v-v -^ 


v'j — v' 


V — V- 


ifiv) - fiv-)) 


2 

fiv'+) - f{v'_) 


+ Re ( Xq 


. v'-v'_ 


V — V- 




Here, we have used Ag + Ag = 1 and 


V — V- 


(±exp (^i (^ipQ - exp (-i(/?g)) = ±7r/2. 


The sign depends on the orientation of the vertices v, v'_v-,v'_^_. But in any case, the expression inside 
arg is purely imaginary. 

Let us write uj = pdz + qdz around Fq or F^, where p,q are functions defined on the vertices of 
Q G R(0) or on the faces incident to u G H(A). Then, by Lemma 12.2.31 and the definition of the 
discrete derivatives 9<}, 9o, 



jj (9o 9 - 9op) Ho 
Fq 

j j [dAq - Bap) Ha 


-2iar(Fg) (d^q - d(^p) 

-2iaT{Fy) {dAq - Bap) 


= j) pdz + j) qdz = j) io, 

OFq BFq BFq 

= ^ pdz + j) qdz = ^ UJ. 

dF^ dFy dFy 


□ 


Note that if w is a discrete one-form of type (), duj = 0 for any face F corresponding to a vertex 
of 0- A discrete one-form uj is said to be closed if do; = 0. Examples for closed discrete one-forms are 
discrete exterior derivatives of complex functions on V (A): 

Proposition 2.3.2. Let f : R(A) — )■ C. Then, ddf = 0. 


Proof. By discrete Stokes’ Theorem 12.3.11 ddf = 0 if d/ = 0 for any discrete elementary cycle P. 
Since df is of type 0, the statement is trivially true if P = Pq for Q G R(0)- So let P = Py for 
V G R(A). Using discrete Stokes’ Theorem 12.3.11 again. 



fjv's) - f{v's-i) 
2 


= 0 . 


□ 

An immediate corollary of the last proposition is the commutativity of discrete differentials, general¬ 
izing the known result for rhombic quad-graphs as provided in [5]. 
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Corollary 2.3.3. Let f : F(0) —^ C. Then, ^o^a/ = d(}d\f. In particular, d\f is discrete holomor- 
phic if f is discrete holomorphic. 

Proof. Due to Proposition 12.3.21 0 = ddf = {d()d\f — d<)d\f) □ 

Remark. Note that even in the generic rhombic case, d\d<)h does not always equal d\d<)h for a 
function h : V{(}) ^ C [5]. Hence, an analog of Proposition 12.3.21 cannot hold for such functions h in 
general. 

Corollary 2.3.4. Let f : P(A) —)• C. Then, f is discrete holomorphic if and only if df = pdz for 
some p : P(0) ^ C. In this case, p is discrete holomorphic. 

Proof. Since all quadrilaterals Q are nondegenerate, the representation of dflgp^ as pdz + qdz is unique 
(see also Lemma [2.3.61 below). Now, df = pAfdz + dAfdz. It follows that / is discrete holomorphic at 
Q if and only if dflgpg = pdz. 

Assuming that df = pdz for some p : V (0) ^ C, ddf = 0 by Proposition 12.3.2] yields d^p = 0. □ 

Let us say that a discrete one-form lo is discrete holomorphic if a; = pdz for some p : P(0) —^ C and 
doj = 0. This notion recurs in the more general setting of discrete Riemann surfaces in the thesis of 
the second author m- By Corollary 12.3.41 df is discrete holomorphic if / is, and by Proposition 12.2.1^ 
any discrete holomorphic one-form w on a simply-connected domain is the discrete exterior derivative 
of a discrete holomorphic function on V (A). 

Due to Chelkak and Smirnov [5], one of the unpleasant facts of all discrete theories of complex 
analysis is that (pointwise) multiplication of discrete holomorphic functions does not yield a discrete 
holomorphic function in general. We can define a product of complex functions on P(A) that is 
defined on V{X) and a product of complex functions on P(A) with functions on R(0) that is defined 
on E{X). In general, the product of two discrete holomorphic functions is not discrete holomorphic 
according to the classical quad-based definition, but it will be discrete holomorphic in the sense that 
a discretization of its exterior derivative is closed and of the form pdz, p : E{X) —)• C, or in the sense 
that it fulfills a discrete Morera’s theorem. 

Corollary 2.3.5. Let f,g : P(A) —C and h : P(0) —>■ C. 

(i) fdg + gdf is a closed discrete one-form. 

(a) If f and h are discrete holomorphic, fhdz is a closed discrete one-form. 

Proof, (i) Let to := fdg + gdf. By Proposition 12.3.21 df and dg are closed. Thus, 

^ a; = f{v) j) dg + g{v) ^ df = 0 
dFy dFy 

for any face E^ corresponding to u E 1^(A). Using Lemma 12.2.31 

2iar(TQ) ^ = 2iar(TQ) j) {fdAgdz -\- fdAgdz + gdAfdz + gdAfdz) 

BFq BFq 

= dAfdAg - dAfdAg + - 5 ac / 9 a / = 0 

for any face Eg corresponding to Q E U(0)- R follows by discrete Stokes’ Theorem 12. 3. II that do; = 0. 
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(ii) By discrete Morera’s Theorem 12. 2. 41 fhdz = 0 for any face F of X, using that / and h are 

discrete holomorphic. Therefore, fhdz is closed. □ 


Remark. In particular, a product f-g'-V (X) —)• C can be defined by integration and f-g is defined up 
to an additive constant. Furthermore, f - h : E{X) —)• C can be defined by “pointwise” multiplication. 
If /, g, h are discrete holomorphic, fdg + gdf = pdz is closed, where p = f ■ d\g + g ■ d\f : E{X) —)■ C, 
and so to say a discrete holomorphic one-form, meaning that f ■ g is discrete holomorphic in this sense. 
Similarly, fhdz is closed, so / • h is discrete holomorphic in the sense that a discrete Morera’s theorem 
holds true. However, / • g and / • h are generally not discrete holomorphic everywhere according to the 
classical quad-based definition of discrete holomorphicity on the dual of a bipartite quad-graph. For 
this, we place any Q (0) as a vertex somewhere in the interior of the corresponding quadrilateral. 

In fact, f ■ g is a complex function on the vertices of X. The medial graph X is not a quad-graph itself, 
but X is the dual of a bipartite quad-graph. More precisely, X is the dual of the bipartite quad-graph 
with vertex set H(A) U H(0), edges connecting points Q E H(0) with all incident vertices v E H(A), 
and faces corresponding to edges of A. Then, f-g does not need to be a discrete holomorphic function 
on the faces of the latter quad-graph. For example, consider /(u) = 0 if u is black and f{v) = lifv 
is white, and a discrete holomorphic g that is not biconstant. Then, the product f ■ g is not discrete 
holomorphic at all Q E H(0) (seen as vertices of the quad-graph described above) where d\g{Q) / 0. 

Furthermore, f ■ h is a complex function on the edges of X, so it is a function on the vertices of the 
medial graph of X. The medial graph of the medial graph of A is usually not a quad-graph, but 
it is the dual of a bipartite quad-graph. Namely, it is the dual of the quad-graph with vertex set 
(H(A) UH(0)) U V{X), edges connecting points v E H(A) or Q E H(0) with the midpoints of all 
incident edges, and each face having an edge of X as a diagonal. Since fhdz is closed, f ■ h is discrete 
holomorphic at vertices of A or <0 by discrete Morera’s Theorem 12.2.41 But there is no need for f ■ h to 
be discrete holomorphic at vertices of X, even for constant h and / defined by f{v) = 0 if u is black 
and f{v) = litv is white. 

In summary, we defined reasonable products / • g and / • h, where f,g:V (A) ^ C and h : V (0) ^ C 
are discrete holomorphic. Somehow missing is a product h ■ h', where h' : H(<C>) — )• C. In the general 
case, we do not know an appropriate product so far. But we want to point out that Chelkak and 
Smirnov defined such a product for so-called spin holomorphic functions h,h' in [^. This product 
satisfies Re (9o {h ■ h'f) = 0. 

2.3.2 Discrete wedge product 

Following Whitney [26], Mercat defined in m a discrete wedge product for discrete one-forms living 
on the edges of A. Then, the discrete exterior derivative defined by a discretization of Stokes’ theorem 
is a derivation for the discrete wedge product. However, a discrete Hodge star cannot be defined on 
A. To circumvent this problem, Mercat used an averaging map to relate discrete one-forms on the 
edges of A with discrete one-forms on the edges of F and F*, i.e., discrete one-forms of type 0- Then, 
he could define a discrete Hodge star; however, the discrete exterior derivative was not a derivation 
for the now heteregoneous discrete wedge product anymore. 

We propose a different interpretation of the discrete wedge product. It the end, we somehow recover 
the definitions Mercat proposed in [IHl HS] [20| , but our derivation is different. Starting with discrete 
one-forms of type 0 that are defined on the edges of X, we obtain a discrete wedge product on the 
faces of X that is of type 0- This definition is different from Whitney’s [26| and has the advantage 
that both a discrete wedge product and a discrete Hodge star can be defined on the same structure. 
In contrast to Mercat’s work, we now can make sense out of the statement that the discrete exterior 
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derivative is a derivation for the discrete wedge product, see Theorem 12. 3. 81 This proposition is of 
crucial importance to deduce discrete integral formulae such as discrete Green’s identities. 

Lemma 2.3.6. Let u he a discrete one-form of type 0- Then, there is a unique representation 
Lo = pdz + qdz with functions p, q : G(0) —^ C. On a quadrilateral Q € 1^(0), P and q are given hy 

P(Q) = ^Q^ + h^. 

Here, e is an edge of X parallel to a black edge ofT, and e* corresponds to a white edge ofT*. 

Proof. Since tv is of type 0, a representation uj = pdz + qdz exists for any quadrilateral Q. In fact, 
given oj, we have to solve a nondegenerate system of two linear equations in the variables p and q. 
Thus, p, q are uniquely defined on G(0)- 

Furthermore, we can find for any quadrilateral Q a function / that is defined on the vertices of Q 
such that 2 f^uj = /(6+) — f(b-) and 2 uj = f{w+) — f{w-), where e is one of the two oriented 
edges of X going from the midpoint of and w± to the midpoint of 6+ and w±, and e* is one of 
the two edges connecting the midpoint of W- and b± with the midpoint of Wj^ and h±. By discrete 
Stokes’ Theorem 12.3.11 we locally get uj = df = pdz + qdz with p = d\f and q = 9a/. Replacing the 
differences of / by discrete integrals of uj yields the desired result. □ 

Definition. Let uj = pdz + qdz and uj' = p'dz + q'dz be two discrete one-forms of type 0, where the 
functions p,p', q, q' ■ R(0) —>■ C are given by Lemma [2.3.61 Then, the discrete wedge product cu A cj' is 
defined as the discrete two-form of type 0 that equals 

{pq - qp') flo 


on faces corresponding to R(<))- 


Note that if one considers dz and dz as discrete one-forms of type 0, 


Ll<) = dz /\ dz. 


Proposition 2.3.7. Let F be a face of X corresponding to Q G R(0), and let e,e* he oriented edges 
of X parallel to the black and white diagonal of Q, respectively, such that Im(e*/e) > 0. Then, 



UJ Auj' = 2 J UJ J uj' — 2 J UJ J ■ 


Proof. Both sides of the equation are bilinear and antisymmetric in uj,uj'. Hence, it suffices to check 
the identity for uj = dz, uj' = dz. On the left hand side, we get JJ^uj Auj' = — 2iar(F). This equals 
the right hand side 


2ee* — 2e*e = 4ilm(ee*) = —i|2e||2e*| sin((/?Q) = —2iai{F). 


□ 

Remark. Since the complex numbers e and e* are just half of the oriented diagonals, we recover the 
definition of the discrete wedge product given by Mercat in [TSl UHl [20] . 
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The discrete exterior derivative is a derivation for the discrete wedge product if one considers functions 
on A and discrete one-forms of type 0: 

Theorem 2.3.8. Let f : F(A) ^ C and let uj he a discrete one-form of type 0- Then, 

d{fio) = df Alo fdio. 

Proof. Let to = pdz + qdz with p,q : V{()) -A C given by Lemma r2.8.6[ If and Fq are faces of X 
corresponding to a vertex v and a face Q of A, 

d{f‘^)\F„ = {fiv) (d^q) (v) - f{v) {d()p) (u)) Qa = fduj\F„, 
d{fuj)\F^ = {q{Q) (9a/) {Q)-p{Q) (9a/) (Q)) = (df A uj)\f^. 

But (df A uj)\f„ = 0 and fduj\FQ = 0, so d(fuj) = df Au -\- fduj. □ 

Remark. In |18] . Mercat proved an analog of Theorem 12.3.81 in a setting where discrete one-forms 
are defined on edges of A. He also considered discrete differential forms as we do, however, the claim 
d(fuj) = df Au + fdoj could not be well defined in his setting. At first, our approach using the medial 
graph allows to make sense out of this statement. It turns out that Theorem 12.3.81 is a very powerful 
theorem leading to discretizations of Green’s identities in Section 12.4.11 and of a Cauchy’s integral 
formula for the discrete derivative of a discrete holomorphic function in Section 12.61 

Above, we defined a discrete wedge product just of two discrete one-forms of type 0- Actually, we 
could define a discrete wedge product of two discrete one-forms of type A in essentially the same way, 
getting a discrete two-form of type A. Then, the analog of Theorem 12.3.81 would be true for this kind 
of discrete wedge product and functions on H(0)- Also the discrete Hodge star of a discrete one-form 
in the next section could be defined not only for those of type <0- However, there exist no analogs of 
Proposition 12.3.71 and 12.3.101 These propositions imply that the discrete wedge product as well as the 
discrete Hodge star of discrete one-forms of type 0 can be defined in a chart-independent way. This 
enables one to consider these objects on discrete Riemann surfaces, see the thesis of the second author 
[12j . There are no such statements if one chooses discrete one-forms of type A. In fact, a discrete 
one-form of type A cannot be canonically defined on a discrete Riemann surface as opposed to discrete 
one-forms of type 0- So since our interest lies in the latter, we do not define a discrete wedge product 
or a discrete Hodge star for discrete one-forms of type A. 


2.3.3 Discrete Hodge star 


Definition. Let / : F(A) <C, h : V(()) —>■ C, w = pdz + qdz a discrete one-form of type <0 with 
complex functions p,q : H(0) —^ C given by Lemma [2.3.61 and Hi,H 2 discrete two-forms of type A 
and 0- The discrete Hodge star is given by 


*/ := -ipf^A] 

2i 


irh —-/iHa; 

2i 


-kuj := —ipdz + iqdz] := —2i 


. Hi 

Hi’ 


aHo :— — 2i 


. H2 
H<^ 


If to and uj' are both discrete one-forms of type 0, we define their discrete scalar product 


(oj,oj') 


iO A , 


F(X) 


whenever the right hand side converges absolutely. Similarly, a discrete scalar product for discrete 
two-forms of the same type is defined. 
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Corollary 2.3.9. (i) = —Id on discrete one-forms of type 0- 

(ii) *2 = Id on complex functions on 1^(A) or 1^(0) and discrete two-forms of type A or 0- 
(Hi) f : F(A) —>■ C is discrete holomorphic if and only if *df = —idf. 

(iv) (/i,/ 2 ) = ffp(x) and = ffjp(x) hi*h 2 for functions /i ,/2 : ^(A) ^ C and func¬ 

tions hi,h 2 '■ V (0) —>• C. 

(v) (•, •) is a Hermitian scalar product on discrete differential forms (of type A or of type ()). 

Proposition 2.3.10. Let uj be a discrete one-form of type (), let Q G 1^(0); and let e,e* be oriented 
edges of X parallel to the black and white diagonal of Q, respectively, such that Im(e*/e) > 0. Then, 

|e-1 JnU)/“’ 

e e e* 

p p* 


Proof. Both sides of any of the two equations are linear and behave the same under complex conju¬ 
gation. Thus, it suffices to check the statement for uj = dz. Hence, it remains to show that 


-ie = cot {(fq) e - 


e’ sm 


(tq) 


e = 


|e|sin((/?Q) 


e - cot (ipq) e* 


Now, both sides of the first equation behave the same under scaling and simultaneous rotation of 
e and e*, the same statement is true for the second equation. Thus, we may assume e = 1 and 
e* = cos {(fq) -t- isin((/?Q). Multiplying both equations by sin(</5Q) gives the equivalent statements 


-isin((/?Q) = cos((/?q) - (cos ((/pq) -|-i sin((/?Q)), 
-isin((^Q)exp(i(/?Q) = 1 - cos ((^q) exp(iv7(3). 


Both equations are true, noting for the second that cos (ipq) — ism{ipq) = exp{—iipq). □ 

Remark. Proposition 12.3.101 shows that our definition of a discrete Hodge star on discrete one-forms 
coincides with Mercat’s definition given in [20]. But on discrete two-forms and complex functions, our 
definition of the discrete Hodge star includes an additional factor of the area of the corresponding face 
of X. 


Proposition 2.3.11. 5 := — -k d-k is the formal adjoint of the discrete exterior derivative d: Let 
f : V{A) — )> C, and let uj be a discrete one-form of type 0 and Ll a discrete two-form of type A. 
Assume that all of them are compactly supported. Then, 

{df,uj) = {f,Suj) and {duj,Q) = {uj,6Ll). 


21 








Proof. By discrete Stokes’ Theorem 12. 3. 11 Theorem I2.3.81 and Corollary 12.3.91 fii). 


0 = JJ d{f-kuj) = JJ dfA-kU} + JJ fd-kUJ = {df,U}) + {f,-kd-kU!), 

F{X) F{X) F{X) 

0 = JJd{-kD,Oj)= JJ {d-kCl) A LV + JJ -kCldoJ = {uj,-kd -k ^) + {doj, Cl). 

F{X) F{X) F{X) 


□ 


2.4 Discrete Laplacian 

The discrete Laplacian and the discrete Dirichlet energy on general quad-graphs were first introduced 
by Mercat in [20]. Later, Skopenkov reintroduced these definitions in [22], taking the same definition 
in a different notation. In our discussion of the discrete Laplacian in Section 12.4.11 we follow the 
classical approach of Mercat (up to sign) and adapt it to our notations. A feature of our notation 
is that we are able to formulate a discrete analog of Green’s first identity and that our formulation 
of discrete Green’s second identity is more intuitive than the previous ones of Mercat, Chelkak and 
Smirnov, and Skopenkov [201 [3122]. 

In Section 12.4.21 we investigate the discrete Dirichlet energy. In particular, we show how uniqueness 
and existence of solutions to the discrete Dirichlet boundary value problem imply surjectivity of the 
discrete differentials and the discrete Laplacian in Theorem 12.4.111 We conclude this section with a 
result concerning the asymptotics of discrete harmonic functions. 

2.4.1 Definition and basic properties 

Definition. The discrete Laplacian on functions / : i4(A) —>■ C, discrete one-forms of type 0, or 
discrete two-forms of type A is defined as the linear operator 

A := —5d — d5 = kdk d + dk dk . 

f is said to be discrete harmonic at u G D(A) if A/(u) = 0. 

The following factorization of the discrete Laplacian in terms of discrete derivatives generalizes the 
corresponding result given by Ghelkak and Smirnov in [3 to general quad-graphs. The local represen¬ 
tation of A/ at u G D(A) is, up to a factor involving the area of the face T), of X corresponding to u, 
the same as Mercat’s in |20j . 

Corollary 2.4.1. Let f : V{A) —)• C. Then, Af = 4dodA/ = Ido^A/. At a vertex v of A, Af{v) 
equals 

E {/W)-/K-i)))- 

Qs^V 

Here, ps := pq„ if v is black, and ps := IIpq„ if v is white. 

In particular, Re(A/) = ARe(/) and Im(A/) = AIm(/). 

Proof. The first statement follows from Af = kdkdf = 2dod\f + 2d(^d\f, noting that d(^d\f = d^dp^f 
by Corollary 12.3.31 
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For the second statement, let us assume without loss of generality that v G F(r). Then, we have to 
show that A/(u) equals 


1 

2ar(F„) 



{f{vs) - fiv)) - COt((y9Qj (/(u') - /(u;_i))^ . 


The structure is similar to the formula of the discrete Hodge star in Proposition 12.3.101 Indeed, if 
denotes an edge of X parallel to the black diagonal vvg and e* an edge parallel to the dual diagonal, 
A/(u) equals 


1 

ar(T;) 




— y 


/ 


ar(F^) 


Qs~V 


\ 


E 


2ar(F^) ^ Vsin((^Q, 


e^l sin(vPQ 
\PQ, 


df-cot{ipQj / df 


ifivs) - fiv)) - coticpqj (/(u') - /(u'_i 



using discrete Stokes’ Theorem 12.3.11 Proposition 12.3.101 and Ipg^l = |e*|/|es|. □ 

Remark. In the case that the diagonals of the quadrilaterals are orthogonal to each other, pq is 
always a positive real number. In this case, the discrete Laplacian splits into two separate discrete 
Laplacians on P and P*. In this case, it is known and actually an immediate consequence of the 
local representation in Corollary 12.4.11 that a discrete maximum principle holds true, i.e., a discrete 
harmonic function can attain its maximum only at the boundary of a region. This is not true for 
general quad-graphs, see for example Skopenkov’s paper [22] . 

Corollary 2.4.2. Let f : P(A) C. 


(i) If f is discrete harmonic, d\f is discrete holomorphic. 

(a) If f is discrete holomorphic, f. Re/, and Im/ are discrete harmonic. 

Proof. By Corollary 12.4.11 A/ = 49 ocIa/ = Ad^dAf. In particular, B^dAf = 0 if A/ = 0, which 
shows (i). Also, / is discrete harmonic if it is discrete holomorphic. Using Re(A/) = ARe(/) and 
Im(A/) = AIm(/), Re(/) and Im(/) are discrete harmonic if / is. □ 

Similar to Proposition 12.2.11 the discrete Laplacian coincides with the smooth one up to order one in 
the general case and up to order two for parallelogram-graphs. This was already shown by Skopenkov 
in [22]. Since this result follows immediately from our previous ones, we give a proof here as well. 

Proposition 2.4.3. Let /c : C ^ C and f its restriction to U(A). 


(i) If fciz) is a polynomial in Re( 2 ;) and Im( 2 ;) of degree at most one, then the smooth and the 
discrete Laplacian coincide on vertices: Acfdv) = A/(u). 

(a) Let all quadrilaterals ofO be parallelograms. If fciz) is a polynomial mRe(z) and Im(z) of degree 
at most two, then the smooth and the discrete Laplacian coincide on vertices: Aqfciv) = A/(u). 
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Proof, (i) It follows from Proposition 12.2.11 (ii) and Corollary 12.4.21 (ii) that constants as well as the 
complex function f(v) = v are discrete harmonic since they are discrete holomorphic. Similarly, 
f{v) = v is discrete harmonic. 

(ii) In the parallelogram case, let us place the vertices of 0 at the centers of the parallelograms. Then, 
f{v) = is discrete harmonic by Proposition 12.2.1] (iiil and Corollary 12.4.21 liil. Looking at real and 
imaginary part separately, A/^ = A/| and A (/ 1 / 2 ) = 0 with fiiy) = Re(u), f 2 {v) = Im(u). Finally, 

A|/P = 4ao4|/P = 49o/i = 4 

with h{Q) = Q, due to Propositions I2.21T] (ivl and 12.2.51 Since any polynomial in Re( 2 :) and Im( 2 :) of 
monomials of degree two is a linear combination of /f — /I; /i + /I i / 1 / 2 , and since we have shown 
that the discrete Laplacian A and the smooth Laplacian Ac coincide on these, we are done. □ 

Remark. The second part of the last proposition generalizes the known result for rhombi given by 
Chelkak and Smirnov in [5]. Note that this is not true for general quadrilaterals even if one assumes 
that the diagonals of quadrilaterals are orthogonal to each other. For this, consider the following 
(finite) bipartite quad-graph of Figure [5J 0 G F is adjacent to the white vertices ±1 and ±z in A, and 
adjacent to the black vertices 2 -|- 2i, — 1 ± i, and 1 — i in F. There are no further vertices. Then, 
A/(0) / 0 for f{v) = v^. Indeed, we would get A/(0) = 0 if we had replaced v = 2 + 2ihyv = l + i 
obtaining a rhombic quad-graph; but |/9Qp/Re(pQ) (/(u) — /(O)) scales by a factor of 2, whereas the 
other nonzero summands in the local representation of A/(0) remain invariant. 



2 + 2i 


Figure 5: A/(0) / 0 for f{v) = 


In the case of general quad-graphs, smooth functions /c : C ^ C, and restrictions / to P(A), 
Skopenkov compared the integral of Ac/c over a square domain R and a sum of A/(u) over black 
vertices of A in R m- Moreover, he showed that for f{v) = |up. 


A/(u) 


2 

ar(F;) 


^ aie&{vv'^_^Qsv'^) 

Qa^V 


when Qs is placed at the intersection point of the middle perpendiculars to the diagonals of the 
corresponding quadrilateral (which equals the intersection point of the diagonals if the quadrilateral is 
a parallelogram). Note that in general, h{Q) = Q is not discrete holomorphic if the vertices Q G P(0) 
are placed on the intersections of these middle perpendiculars. 


For a finite subset Oo C 0 and two functions f,g : P(Ao) —>■ C, we denote by 

{f,9)Oo ■=-Yi // 

F{Xo) 
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the discrete scalar product of / and g restricted to Oo- Similarly, the restriction of the discrete scalar 
product of two discrete one-forms of type 0 on Oo is defined. 

In the rhombic setup, discrete versions of Green’s second identity were already stated by Mercat 
[18j . whose integrals were not well defined separately, and Chelkak and Smirnov [5], whose boundary 
integral was an explicit sum involving boundary angles. Skopenkov formulated a discrete Green’s 
second identity with a vanishing boundary term m- We are able to provide a discrete Green’s first 
identity for the first time, prove it completely analogonsly to the smooth setting, and dednce discrete 
Green’s second identity out of it. 

Theorem 2.4.4. Let Oo 0 be finite, and let f,g : y(Ao) —>■ C. 

(i) {f,^g)0o + idf,d9)0o= f f-*^dg. 

axo 

("Aj (A/,5r)oo - (/, A5r)oo = f {f-kdg - g-kdf). 

aXo 


Proof, (i) Since the discrete exterior derivative is a derivation for the discrete wedge product by 
Theorem 12.3.81 d{f * dg) = df A -kdg + / * {kd * dg). Now, discrete Stokes’ Theorem 12.3.11 yields the 
desired result. 

(ii) Just apply twice discrete Green’s first identity. □ 

Discrete Weyl’s lemma is a direct consequence of discrete Green’s second identity, Theorem 12.4.41 (ii). 
A version for rhombic quad-graphs was given by Mercat in [18], proven by an explicit calculation. 

Corollary 2.4.5. / : 14(A) —>■ C is discrete harmonic if and only if {f,Ag) = 0 for every compactly 
supported g : V (A) ^ C. 

Skopenkov introduced the notion of discrete harmonic conjugates in |22]. We recover his definitions in 
onr notation, observing that his discrete gradient corresponds to our discrete exterior derivative and 
his counterclockwise rotation by 7r/2 corresponds to our discrete Hodge star. 

Definition. Let Oo ^ 0 and / be a real (discrete harmonic) function on 14(Ao). A real discrete 
harmonic function / on 14(Ao) is said to be a discrete harmonic conjugate of f ii f + if is discrete 
holomorphic. 

Note that the existence of a real fnnction / such that f + if is discrete holomorphic requires already 
that / is discrete harmonic due to Corollary 12.4.21 (ii). 

Lemma 2.4.6. Let Oo ^ 0 a,nd f be a real discrete harmonic function on 14(Ao). 

(i) The discrete harmonic conjugate f is unique up to two additive real constants on Lq and Lg. 

(ii) If Oo is simply-connected, a discrete harmonic conjugate f exists. 

Proof, (i) If /i and /2 are two real discrete harmonic conjugates, their difference /i — /2 is real and 
discrete holomorphic. So by Proposition 12.2.2] fiil. it is biconstant. 

(ii) Since / is harmonic, d k df = 0, i.e., kdf is closed and of type 0- Moreover, reality of / implies 
kdf = —id^fdz -p id\fdz = 2 Im {d\fdz). So in the same manner as in the proof of Proposition 12.2.81 
kdf can be integrated to a real function / on 14(Ao). Finally, f + if is discrete holomorphic by 
Corollary 12.3.41 since d{f + if) = df + ikdf = 2Re {d\fdz) -|- 2i Im [d/^fdz) = 2di^fdz. □ 
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Note that in the case of quadrilaterals with orthogonal diagonals, such that A splits into two discrete 
Laplacians on T and F*, it follows that a discrete harmonic conjugate of a discrete harmonic function 
on l^(r) can be defined on F(r*) and vice versa, as was already noted by Chelkak and Smirnov in [5]. 

Corollary 2.4.7. Let Oo C 0 and f : l^(Ao) —>■ C 6e discrete holomorphic. Then, Im(/) is uniquely 
determined by Re(/) up to two additive real constants on Fq and Fq. 


2.4.2 Discrete Dirichlet energy 


We follow the classical approach of discretizing the Dirichlet energy introduced by Mercat in [20] . Note 
that Skopenkov’s definition in [22] is exactly the same. In particular, Skopenkov’s results, including an 
approximation property of the Laplacian, convergence of the discrete Dirichlet energy to the smooth 
Dirichlet energy for nondegenerate uniform sequences of quad-graphs, and further theorems for quad- 
graphs with orthogonal diagonals apply as well in our setting. We refer to his work [22] for details on 
these results. 


Definition. Let Oo ^ 0- For / ; D(Ao) —C, we define the discrete Dirichlet energy of / on Oo as 
^Oo(/) := (4f,4f)Oo £ [0,oo]. 

If Oo is finite, the discrete Dirichlet boundary value problem asks for a real discrete harmonic function 
/ on D(Ao) such that / agrees with a preassigned real function /o on the boundary D(9Ao). 


Proposition 2.4.8. 


EoM) 


E 

QGV(Oo) 


1 

2Re(pQ) 


{\PQ? !/(&+) - + \f{w+) - /(u;_)|^) 


+ E (^(/(6+) -/(6_)) (/(u;+) -/(u;-))) . 

QeV(Oo) 


Proof. Let us denote the right hand side by Then, Sog(/) = E(^^(Re{f)) -|- Sog(Im(/)) and 

E'^^{f) = £'^^(Re(/)) -|- ii'^^(Im(/)). Therefore, we can restrict to real functions /. Furthermore, it 
suffices to check the identity for just a singular quadrilateral Q, because both F^0(,(/) and E'^^{f) are 
sums over Q € R(Oo)- 

Also, both expressions depend only on df, so we may assume f{b-) = f{w-) = 0. In addition, both 
terms do not change if Q is rotated or scaled. Thus, we can assume 6+ — 6_ = 1 as well, such that 
t(;+ — w- = ipQ. Then, 


EqU) = / / df A-kdf = 4area((3)|aA/| = 2Re(pQ)|aA/|^ 


Fq 

=2Re{pQ) 


>^QfiK) + —fiw+) 

I'PQ 


=2|pQ|sin((/?Q) 




+ 




4sin2 ((^ Q ) As\v?{pq)\pq\^ 


2Re(pQ) 


r{K) + 


2Re(/9Q) 


/ {w+) + 2f{b+)f{w+) 


+ 2/(6+)/(u;+)Im 

Im { pq ) 




2Re(/3Q) 


fQj^ 

= EqU)- 


Here, we used that 2area(Q) = |t(;+ - r(;_| sin((/?Q) = \pq\ sin((/?Q) = Re(pQ), 1 /|Aq| = 2sin((/9Q), and 
Aq = —/9Q/(4area((5)) (see the proof of Lemma [2. 2. 3h . □ 
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The same formula of i? 0 p(/) was given by Mercat in |20j . 

In the case of rhombic quad-graphs, Duffin proved in [9] that the discrete Dirichlet boundary value 
problem has a unique solution. The same argument applies for general quad-graphs with the discrete 
Dirichlet energy defined here. Using a different notation, Skopenkov proved existence and uniqueness 
of solutions of the discrete Dirichlet boundary value problem as well [22]. 

Lemma 2.4.9. Let Oo C 0 finite and fo : U((9Ao) —>■ IR. We consider the vector space of real 
functions f : U(Ao) — IR that agree with /o on the boundary. 

Then, is a strictly convex nonnegative quadratic functional in terms of the interior values f{v). 
Furthermore, 

-^(/) = 2ar(F,)A/M 

for any v E U(Ao)\U(5Ao). In particular, the discrete Dirichlet boundary value problem is uniquely 
solvable. 

Proof. By construction, £<}(, is a quadratic form in the vector space of real functions / : V (A) —?■ IR. In 
particular, it is convex, nonnegative, and quadratic in terms of the values f{v). Thus, global minima 
exist. To prove strict convexity, it suffices to check that the minimum is unique. 

For an interior vertex vq E U(Ao)\U((9Ao), let 4){v) := d^^vo be the Kronecker delta function on U(Ao). 
Then, 

= ^Eofif + t(j))\t=o = 2((i/,#)oo = -2{A f,(j)) = -2ai{Fy^)Af{vo) 

due to Proposition 12.3.111 It follows that minima of FIoq are discrete harmonic in the interior of Oo- 
The difference of two minima is a discrete harmonic function vanishing on the boundary, so its energy 
has to be zero. But only biconstant functions have zero energy. Thus, the difference has to vanish 
everywhere, i.e., minima are unique. □ 

In the following, we apply Lemma [2.4.91 to show that dA,d\,do,d(), A are surjective operators. This 
implies immediately the existence of discrete Green’s functions and discrete Cauchy’s kernels, as we 
will see in Sections 12.51 and 12.61 

Lemma 2.4.10. Let ()o T 0 be finite and homeomorphic to a disk. Then, the discrete derivatives 
d\,d\,do,d() and the discrete Laplacian A are surjective operators. That means, given any complex 
functions ho on U(Oo) o.nd fo on U(Ao\cIAo), there exist functions hQ,hQ on U(Oo) o.nd functions 
fd, fd, /a on U(Ao) such that d^hQ = B^^hg = A/a = fo and Sa/s = dAfg = ho. If fo is real-valued, 
/a can be chosen real-valued as well. 

Proof. Let V, E, and E denote the number of vertices, edges, and faces of Aq, respectively. Denote by 
BV and BE the number of vertices and edges of SAq, respectively. By assumption, cIAq is a simple 
polygon and BV = BE. 

By Lemma 12.4.91 the space of real discrete harmonic functions on U(Ao) has dimension BV. Clearly, 
real and imaginary part of a discrete harmonic function are itself discrete harmonic. Therefore, the 
complex dimension of the space of complex discrete harmonic functions, i.e., of the kernel of A, is 
BV as well. Thus, we have shown that A : —>■ is a surjective linear operator with 

K E {R,C}. 

Now, A = 4:dodA = 49of?A by Corollary 12.4.II so —)■ ^re surjective as well. 

The kernel of B() consists of all discrete holomorphic functions on U(Oo)- By Proposition 12.2.^ any 
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such function has a discrete primitive, i.e., the kernel is contained in the image of d\. Using the 
surjectivity of A, it follows that 9 a : is surjective. The same is true for d\. □ 

Theorem 2.4.11. The discrete derivatives d\,d\,d(),d() and the discrete Laplacian A (defined on 
complex or real functions) are surjective operators on the vector space of functions on U(A) or U(0)- 


Proof. Let Oo C 0i C O 2 C ... C 0 be a sequence of simply-connected finite domains such that 
U^o = 0- By Afc we denote the subgraph of A whose vertices and edges are the vertices and edges 
of quadrilaterals in Ofc- 


Let us first prove that any h : U(0) —^ C has a preimage under the discrete derivatives 9 a, 9a. By 
Lemma [2.4.101 the affine space A® of all complex functions on V{Ak) that are mapped to ^|y( 0 j,) 


by 9a (or 9a) is nonempty. Let 
V{Aj) C V{Ak). Clearly, 


A, 


denote the affine space of restrictions of these functions to 


D a)' 


( 0 ) 


Ao 


D A2 


(0) 


D 


Ao 


Since all affine spaces are finite-dimensional and nonempty, this chain becomes stationary at some 
point, giving a function /o on U(Ao) mapped to by d\ (or d\) that can be extended to a 

function in A^^^ for any k. 

Inductively, assume that fj : V{Aj) —>■ C is mapped to by d\ (or 9a) and that fj can be 

extended to a function in A^^^ for all k > j. Let k > j + 1, he the affine space of all complex 

functions on U(Afc) that are mapped to /i|y( 0 j.) by d\ (or 9a) and whose restriction to V{Aj) is equal 
to fj. By assumption, all these spaces are nonempty. In the same way as above, there is a function 
/j+i extending fj to U(Aj+i) that is mapped to by d\ (or 9a) and that can be extended to 

a function in for sd\ k > j + 1. 

For v G U(Ai.), define fiv) := fkiv). f is a well defined complex function on U(A) with d\f = h (or 
9a/ = h). Hence, 9a, 9a : ^ are surjective. 

Replacing U(C>fe) by V{Ak\dAk), we obtain with the same arguments that A is surjective, regardless 
whether A is defined on real or complex functions. 

Finally, 9o,9^ : —>• are surjective due to A = 49o9a = 49o9a by Corollary 1 2.4.1 1 □ 


In the case of rhombic quad-graphs with bounded interior angles, Kenyon proved the existence of 
a discrete Green’s function and a discrete Cauchy’s kernel with asymptotic behaviors similar to the 
classical setting m- But in the general case, it seems to be practically impossible to speak about any 
asymptotic behavior of certain discrete functions. For this reason, we will not require any asymptotic 
behavior of discrete Green’s functions and discrete Cauchy’s kernels in Sections 12.51 and 12.61 

Still, one can expect results concerning the asymptotics of special discrete functions if the interior 
angles and the side lengths of the quadrilaterals are bounded, meaning that the quadrilaterals do not 
degenerate at infinity. And indeed, on such quad-graphs we can show that any discrete harmonic 
function whose difference functions on U(r) and U(F*) have asymptotics o(u“^/^) is biconstant. In 
the rhombic setting, Chelkak and Smirnov showed that a discrete Liouville’s theorem holds true, i.e., 
any bounded discrete harmonic function on V (A) vanishes [5] . 

Theorem 2.4.12. Assume that there exist constants oq > 0 and 00 > Ei > Eq > 0 such that a > ao 
and El > e > Eq for all interior angles a and side lengths e of quadrilaterals Q & V (0) • 
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_ 2 /2 

U f ■ ^ C is discrete harmonic and /(w+) — f{v-) = o{v^ ' ) for any two adjacent v± G l^(r) 

or v± € / is biconstant. 

Proof. Without loss of generality, we can restrict to real functions /. Assume that / is not biconstant. 
Then, df A -kdf is nonzero somewhere on a face F of X. In particular, the discrete Dirichlet energy of 
/ is bounded away from zero if the domain contains F. Now, the idea of proof is to show that if the 
domain is large enough, the function being zero in the interior and equal to / on the boundary has a 
smaller discrete Dirichlet energy than /, contradicting Lemma 12.4.91 

Let us first bound the intersection angles and the lengths of diagonals of the quadrilaterals. Take 
Q G 1^(0)- Then, there are two opposite interior angles that are less than tt, say a± at vertices h±. 
Since all interior angles are bounded by oq from below, one of a± is less than or equal to n — oq, say 
ao < Oi- ^ — «o- 

By triangle inequality, |6+ — 6_|, |t(;+ — W-\ < 2Ei. Twice the area of Q equals 

|r(;_ — 6_||r(;+ — 6_| sin(a_) + |tc_ — 6+||r(;+ — 6+| sin(a+) > Eq sin(Q;o). 

On the other hand, twice the area of Q is equal to |6+ — 6_||r(;+ — rc-l sin((/9Q). Combining the two 
inequalities yields |6+ — b-\,\w+ — w-\ > Eq := Eq sm{ao)/{2Ei). Furthermore, we have the lower 
estimate sin((/;Q) > Eq sm{aQ)/(iEf) = Eq/{2Ei). Thus, we can bound 

PQ = exp (i (^<fQ - I) ) = (sm(v9Q) - i cos((^q)) 

by \pq\ < ^ and Re(pQ) > • 

For some r > 0, denote by Bo(0, r) C V (0) the set of quadrilaterals that have a nonempty intersection 
with the open ball around 0 and radius r. Let R > 2Ei, and consider the ball i?o(0, i?) C D(0)- 

Since edge lengths are bounded by Ei, all faces in the boundary of i?o(0, R) are contained in the set 
B{0,R + 2Ei)\B(0, R — 2Ei) C C. The area of the latter is SttREi. Any quadrilateral has area at 
least Fig sin(ao)/2, so at most 167ri?Fii/(£'Q sin(Q;o)) quadrilaterals are in the boundary of B(^{0,R). 

Consider the real function fn defined on the vertices of all quadrilaterals in B<){0,R) that is equal to 
/ at the boundary and equal to 0 in the interior of i?o(0, R). When computing the discrete Dirichlet 
energy, only boundary faces can give nonzero contributions. If we look at the formula of the discrete 
Dirichlet energy in Proposition 12.4.81 and use that /(u+) — /(u_) = at the boundary, we see 

that any contribution of a boundary face has asymptotics o{R~^). For this, we use that |Re(/3Q))| 
is bounded from below by a constant and |Im(/?Q)| < \pq\ < 2E\IE'q. Using that there are only 
0{Rf) faces in the boundary, the discrete Dirichlet energy F'_B^(o,R)(/i?)) considered as a function of R, 
behaves as o(l). So if R is large enough, 

EB0{o,R)ifR) < jj df A-kdf < Rb^(o,r)(/), 

F 

contradicting that / minimizes the discrete Dirichlet energy by Lemma 12.4.91 □ 

2.5 Discrete Green’s functions 

In this section, we discuss discrete Green’s functions. 
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Definition. Let vq G A real function G{-;vo) on 1^(A) is a (free) discrete Green’s function for 

Vo if 


G(uo;uo) = 0 and AG(u;uo) = 


1 


-5^0 all V € 1^(A). 


2ar(F^g) 

Corollary 2.5.1. A discrete Green’s function exists for any vq G 1^(A). 


Proof. By Theorem 12.4.111 there exists G : 1^(A) —>■ IR such that AG{v) = (4ar(F„g)). Since 

constant functions are discrete harmonic, we can adjust G to get G{vo) = 0. □ 

Remark. In contrast to the smooth setting or the rhombic case investigated by Kenyon m and 
Chelkak and Smirnov [5], we do not require any asymptotic behavior of the discrete Green’s function. 
But when considering planar parallelogram-graphs with bounded interior angles and bounded ratio of 
side lengths in Section [3.21 of Section O we will prove the existence of a discrete Green’s function with 
asymptotics generalizing the corresponding result for rhombic quad-graphs. 


Our notion of discrete Greens’ functions in a discrete domain and the proof of their existence follow 
the presentation of Ghelkak and Smirnov in [5]. 

Definition. Let Oo C 0 be finite. For vq G K(Ao\9Ao), a real function Gao(s^^o) on K(Ao) is a 
discrete Green’s function in Aq for vq if 

Gao{v; Vq) = 0 for all v G K(9Ao) 

and AGao(v;vo) = i—7Tr-r<^vvo for all v G K(Ao\(9Ao). 

4:31 [r ) 

Corollary 2.5.2. Let Oo C 0 be finite and vq G K(Ao\9Ao). Then, there exists a unique discrete 
Green’s function in Aq for vq . 

Proof. Let Gq : V (A) ^ C be the free discrete Green’s function for vq given by Gorollary 12.5.11 By 
Lemma [2.4.91 there exists a real discrete harmonic function Gi on V (Aq) such that Gi and Gq coincide 
on K(9Ao). Thus, Gao{-',vo) := Gq — Gi is a discrete Green’s function in Aq for vq. Since the difference 
of two discrete Green’s functions in Aq for vq is discrete harmonic on K(Ao) and equals zero on the 
boundary l/(9Ao), it has to be identically zero by Lemma 12.4.91 □ 


2.6 Discrete Cauchy’s integral formulae 

In this section, we first formulate discretizations of the standard Gauchy’s integral formula, both for 
discrete holomorphic functions on V{A) and K(0)- Later, we give with Theorem 12.6.31 a discrete 
formulation of Gauchy’s integral formula for the derivative of a holomorphic function. We conclude 
this part with Section [2.6.11 where we relate our formulation of the discrete Cauchy’s integral formula 
applied to a discrete holomorphic function on V (A) with the notation used by Chelkak and Smirnov 
in [5]. 

Definition. Discrete Cauchy’s kernels with respect to Qq G K(0) o,nd vq G K(A) are functions 
Kqq : V{k) —>■ C and : K(0) — >■ C, respectively, that satisfy for all Q G K(0), v G V{k): 

TT - TT 

dAKqfiQ) = 6 qqo and d(^K^fiv) = 

oxyrg) ar(i<^j 
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Remark. As for the discrete Green’s functions in Section 12.51 we do not require any asymptotic 
behavior of the discrete Cauchy’s kernels. However, if interior angles and side lengths of quadrilaterals 
are bounded, it follows from Theorem 12.4.121 that any discrete Cauchy’s kernel with respect to a 
vertex of 0 with asymptotics is necessarily unique. In Section [3.31 we will construct discrete 

Cauchy’s kernels with asymptotics similar to the smooth setting, generalizing Kenyon’s result m on 
rhombic quad-graphs to parallelogram-graphs. 

The existence of discrete Cauchy’s kernels follows from Theorem l2.4.11l 

Corollary 2.6.1. Let Qo G 1^(0) and vq G K(A) be arbitrary. Then, discrete Cauchy’s kernels with 
respect to Qo and vq exist. 

Theorem 2.6.2. Let f and h be discrete holomorphic functions on V{A) and K(0); respectively. Let 
Vq G K(A) and Qo G K(0), and let Ky^ : V{()) —>• C and Kqq : K(A) C be discrete Cauchy’s kernels 
with respect to vo and Qo, respectively. 

Then, for any discrete contours Cy^ and Cg^ on X surrounding vo and Qo, respectively, once in 
counterclockwise order, discrete Cauchy’s integral formulae hold: 

f{vo) = ^ ^ fKy^dz, 

CvQ 

HQo) = ^ y ^^Qodz- 

CQo 

Proof. Let Py and Pq be discrete elementary cycles, v being a vertex and Q a face of A. By 
Lemma 12.2.31 and the definition of d(), we get: 


1 

2TTi 


f 


Pv 


fKy^dz = -av{Fy)f{v)d(;,Kyg{v) 
7T 


dvvof{v) ) 


7^ / fKy^dz = -ai{FQ)dKf{Q)Ky^{Q) = 0. 

Im J vr 

Pq 


By definition, the discrete contour Cy^ bounds a topological disk, and we can decompose the integration 
along Cyg into a couple of integrations along discrete elementary cycles Py and Pq as above. Summing 
up, only the contribution of Py^ is nonvanishing, and we get the desired result. The second formula is 
shown in an analog fashion. □ 


Remark. In the case of rhombic quad-graphs, Mercat formulated a discrete Cauchy’s integral formula 
for the average of a discrete holomorphic function on K(A) along an edge of A. In [5], Chelkak and 
Smirnov provided a discrete Cauchy’s integral formula for discrete holomorphic functions on K(0) 
using an integration along cycles on T and T*, see Section 12.6.11 

Theorem 2.6.3. Let f : K(A) <C be discrete holomorphic, Qo G K(0), and let Kqq : K(A) —>■ C 6e 
a discrete Cauchy’s kernel with respect to Qq. 

Then, for any discrete contour Cqq in X surrounding Qq once in counterclockwise order that does not 
contain any edge inside Qo (see Figure W, the discrete Cauchy’s integral formula is true: 

dAfiQo) = y fdAKqgdz. 

Cqo 


31 










Proof. Let D be the discrete domain in X bounded by Cq^. Since no edge of Cq^ passes through Qq, 
the discrete one-form df^Kq^dz vanishes on Cq^. Therefore, 


y" fdAKq.dz = j> fdKq, = jj d{fdKq,) = JJ df A dKq, 

Cqo Cq^ d d 

due to discrete Stokes’ Theorem 12.3.11 Theorem 12.3.81 and Proposition 12.3.21 Now, / is discrete 
holomorphic, so df AdKq^ = 9a/9a-K’qo^ 0- I^^t d\Kq^ vanishes on all vertices of 0 but Qq. Finally, 

^ ~^i jj ^AfdAKqo^o = dAf{Qo). 

CQo ^Qo 


□ 


0--f--0-# 



Figure 6: Discrete contour as in Theorem 12.6.31 


Remark. In general, there exists no analog of Theorem 12.6.31 for the discrete derivative of a discrete 
holomorphic function on 1^(0), because the discrete derivative itself does not need to be discrete 
holomorphic. However, in the special case of integer lattices, any discrete derivate of a discrete 
holomorphic function is itself discrete holomorphic. In Section 13.41 we will obtain discrete analogs of 
Cauchy’s integral formulae for higher derivatives of discrete holomorphic functions on P(A) or H(<C>)- 

2.6.1 A different notation 

Let VF be a cycle on the edges of F*, having (counterclockwise ordered) white vertices 11 ) 0 , 101 ,, Wm, 
Wm = vjQ. Then, any edge connecting two consecutive vertices Wk,Wk+i forms the diagonal of a 
face Q{wk,Wk+i) G 1^(0)- We denote the set of such faces together with the induced orientation 
of their white diagonals by W(). For Q G VFo, we denote its white vertices by W-{Q),w+{Q) such 
that the corresponding oriented diagonal goes from W-{Q) to Wj^{Q). Its black vertices are denoted 
by h{Q),h'{Q) in such a way that W-{Q),b{Q),w+{Q),b'(Q) appear in counterclockwise order. The 
reason why we do not choose our previous notation of Figure [His that black and white vertices now 
play a different role that shall be indicated by the notation. 

Now, we construct a cycle B on the edges of F having (counterclockwise ordered) black vertices 
bo,bi,... ,bn, bn = bo, in the following way. We start with bo := b{Q {wo,wi)). In the star of the 
vertex tci, there are two simple paths on F connecting bo and h{Q {wi,W 2 )), and we choose the path 
that does not go through Q{wo,wi). Note that it may happen that b {Q {wi,W 2 )) = bo', in this case, 
we do not add any vertices to B. Also, W 2 = wo is possible, which causes adding the nondirect path 
connecting bo and b{Q { 101 , 102 )) = b' {Q {wo,wi)). 
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Continuing this procedure till we have connected b {Q {wm-i, Wm)) with Bq, we end up with a closed 
path B on F. Without loss of generality, any two consecutive vertices in B are different. As above, any 
edge connecting two consecutive vertices bk,bk+i forms the diagonal of a face Q{bk,bk+i) G 1^(0)- We 
denote the set of such faces together with the induced orientation of their black diagonals by i?o. For 
Q G i?o, we denote its black vertices by b-{Q),b^{Q) such that the corresponding oriented diagonal 
goes from b-{Q) to 6+((5). Finally, its white vertices are denoted by w{Q),w'{Q) in such a way that 
b^{Q),w'{Q),b+{Q),w{Q) appear in counterclockwise order. 


Definition. Let W and B be as above and h a function dehned on W(^ U B(^. We define the discrete 
integrals along W and B by 


m—1 


f h{Q)dz ■='^ f [Q {wk,Wk+i)) (wk+i - Wk ), 

k=0 

~ m—1 

(f) h{Q)dz := ^ / (Q {wk,Wk+i)) {wk+i - Wk)] 
JW k=Q 

» n—1 

j) h{Q)dz := ^ / (Q {bk, bk+i)) (bk+i - bk ), 

k=0 

~ n—1 

(I h{Q)dz :='^f {Q {bk, bk+i)) {bk+i - bk). 

B u—c\ 


In between the closed paths B and W^ there is a cycle P on the medial graph X that comprises exactly 
all edges [Q, n] with Q G kF<) and v & B incident to Q and all edges [Q, v] with Q G So and n G VF 
incident to Q. The orientation of [Q,v] is induced by the orientation of the corresponding parallel 
white or black diagonal. Figure [7] gives an example for this construction. 



Figure 7: Cycles W on F*, B on F, and closed path P on A in between 


Remark. Note that a discrete contour P on A induces a white cycle W and a black cycle B in such 
a way that W, P, and B are related as above. 


Lemma 2.6.4. Let W, B, and P be as above. Let f and h he functions defined on the vertices ofW 
and B and on VFo C B(^, respectively. Then, 


/ 

Jw 


f{b{Q))h{Q)dz + 


<j> f{w{Q))h{Q)dz = 2 ^ fhdz. 


33 




Proof. Any edge e = [Q,b{Q)] {Q E W(f) or [Q,w{Q)] {Q E -Bo) of P corresponds to either an edge 
W-{Q)w+{Q) of W or to an edge b-{Q)b^{Q) of B, respectively, and vice versa. Since this edge of 
X is just half as long as the parallel edge of F or F*, 2 fh = f {b (Q)) h{Q){w+ — W-) in the first 
and 2j^fh = / (rc (Q))/i(Q)(6_|_ — 6_) in the second case. Therefore, the discrete integral along P 
decomposes into one along B and one along W. □ 


Note that the construction of B and Lemma 12.6.41 are also valid if W consists of a single point or of 
only two edges (being the same, but traversed in both directions). In any case, P will be a discrete 
contour. 

The discrete Cauchy’s integral formula in the notation of Chelkak and Smirnov in [5] reads as 


TTih{Qo) = (f h{Q)KQ^{b{Q))dz + (f h{Q)KQ^{w{Q))dz 
Jw JB 

if Qo E C(0) is surrounded once by VF, h is discrete holomorphic on F(0), and Kqq : F(A) — >■ C is a 
discrete Cauchy’s kernel with respect to Qq. Lemma 12.6.41 directly relates this formulation to ours in 
Theorem 12.6.21 

Finally, we conclude this section with a proposition relying on the decomposition of a discrete contour 
into black and white cycles. In Corollary 12.8.51 (i), we have already seen that fdg + gdf is closed for 
functions f,g : F(A) —>■ C. Actually, a slightly stronger statement is true: 

Proposition 2.6.5. Let W and B be as above, and let f,g : F(A) —)• C. Then, 


(f f{b{Q)) {d\g{Q)dz + dAg{Q)dz) + (f g{w{Q)) {dAf{Q)dz + dAf{Q)dz) = 0. 

Jw Jb 

Proof. We first rewrite the discrete integral along W. By definition, we have dg = dAgdz + dAgdz. 
Using discrete Stokes’ Theorem 12.3.11 

/ f{b{Q)){dAg{Q)dz + dAg{Q)dz) = V f{h{Q)){g{w+{Q)) - g{w-{Q))). 

Q€W, 

We can change the summation along the path W into a summation along the path B by 

X] fiKQ)) {a {w+ {Q)) - g {w- (Q))) = - g{w (Q)) if (6+ (Q)) - / (6_ (g))). 

QeW(, QeB^:, 


In an analog way as for W, we can rewrite the discrete integral along B as 



g(w(Q)) (dAf(Q)dz + dAf{Q)dz) = Y, 9{wm (/ (6+ {Q)) - f (6- (Q))) ■ 

Q(zBo 


In summary, 


/ 

Jw 


f{b{Q)) {dAg{Q)dz + dAg{Q)dz) + 



g{w{Q)) [dAf{Q)dz + dAf{Q)dz) 


= 0 . 


□ 
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3 Discrete complex analysis on planar parallelogram-graphs 


This part deals with discrete complex analysis on planar parallelogram-graphs A that are assumed to 
be locally finite. The vertices Q of the dual 0 are placed at the centers of the parallelograms Q. Note 
that any planar parallelogram-graph is strongly regular and bipartite. As in Section [2l the induced 
graphs on black and white vertices are denoted by T and T*, respectively. 

In Propositions 12.2Tl 12.2.51 and 12.4.31 we have already seen that discrete complex analysis on parallel¬ 
ogram-graphs is closer to the classical theory than on general quad-graphs. For example, f{v) = 
is a discrete holomorphic function on 14(A) and d\f{Q) = 2Q; h{Q) = Q is a discrete holomorphic 
function on 14(0) and d()h = 1; and the discrete Laplacian A approximates the smooth one correctly 
up to order two. 

In order to concentrate on the calculation of the asymptotics of a certain discrete Green’s function 
and discrete Cauchy’s kernels, we postpone the discussion of some necessary combinatorial and ge¬ 
ometric results on planar parallelogram-graphs to the appendix. Our setup is closely related to the 
quasicrystallic parallelogram-graphs discussed by the first author, Mercat, and Suris [T]. In their case, 
the quad-graph can be lifted to and a discrete exponential defined on can be restricted to the 
quad-graph. We adapt their construction to our setting in Section 13.11 The discrete exponential is 
the basic building block for the construction of a discrete Green’s function in Section [3.21 and discrete 
Cauchy’s kernels in Section [331 These constructions base on the ideas of Kenyon in the case of rhombic 
quad-graphs m- 

The corresponding functions can be defined for general planar parallelogram-graphs, but we need more 
regularity of the graph to calculate their asymptotics. The two conditions we use are that all interior 
angles of the parallelograms are bounded (the same condition was used in the presentation of Chelkak 
and Smirnov in [5]) and that the ratio of side lengths of the parallelograms is bounded as well. For 
rhombic quad-graphs, the second condition is trivially fulfilled; for quasicrystallic graphs, there are 
only a finite number of interior angles. Note that instead of using boundedness of the ratio of side 
lengths of the parallelograms, we can assume that the side lengths themselves are bounded. This seems 
to be a stronger condition at first, but actually, both conditions are equivalent, see Proposition IA.3I in 
the appendix. 

We conclude this part by a discussion of integer lattices in Section 13.41 On these graphs, discrete 
holomorphic functions can be discretely differentiated infinitely many times, and for all higher order 
discrete derivatives, discrete Gauchy’s integral formulae with the right asymptotics hold true. 

During this section, we use the following shorthand notation. 

Definition. Let v^v’ G 14(A) and Q.,Q' G 14(0)- 

(i) Choose any directed path of edges ei,..., e„ on A going from v' to v. Define 

n 

J{v,v') 

(ii) Choose any directed path from v to Q that begins with a directed path on A with edges ei,..., e„ 
to a vertex vq of the parallelogram Q and that ends with two half-edges dij2,d2l2, where di,d 2 
are emanating from vq. Define 

n 

-J{v,Q) = J{Q,v) 

i=i 
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Moreover, let t{v,Q) = t{Q,v) := l/{did2) if vq,v are both in l^(r) or both in y(r*) and 
t{v,Q) = t{Q,v) := —l/{did2) otherwise. 

(iii) Choose any directed path from Q' to Q consisting of half-edges ei/2, 62/2 connecting the center 
of the parallelogram Q' to one of its vertices vq/, a directed path on A with edges 63 ,... ,e„ 
going from vq/ to a vertex vq of the parallelogram Q, and two half-edges di/2, ^ 2/2 emanating 
from Vq. Define 

J{Q,Q') ■= -e^ ^ + -63 ^ ^ ^ 

i=3 

Furthermore, let t{Q,Q') := 1 /( 6162 ^ 1 ^ 2 ) if vq,vqi are both in D(r) or both in l/(r*) and 
t{Q,Q') := — 1 /( 6162 ^ 1 ^ 2 ) otherwise. 

It is easy to see that these definitions do not depend on the choice of paths. 

Remark. In the case that all parallelograms are rhombi of side length one, we have J(x, x') = x — x'. 


3.1 Discrete exponential function 


Definition. Let vq G i^(A). Then, the discrete exponentials e(-, •; uq), exp(-, •; uq) ^ C x l^(^) C are 
defined by 


e{X,vo]Vo) = 1 = exp(A,i;o;uo), 
e(A, v'] Vq) _ X + (v' - v) 
e{X,v,vo) X-{v'-v)' 
exp(A,u^;uo) ^ 1 ^{v' - v) 

exp{X,v;vo) 

for all vertices v,v' € i^(A) adjacent to each other and all A G C. 

For a face Qo G 1^(0) with incident vertices v-,v'_,v+,v'j^ in counterclockwise order and vq G i^(A), 
we define the discrete exponentials e(-, •; Qo) : C x i^(A) —>■ C and exp(-, •; vq) : C x 1^(0) —>■ C as 


e(A,u;Qo) 

exp(A,(5o;?^o) 


_e(A,u;u±)_ 

(A - (u± - u^)) (A - (u± - u^)) ’ 
exp(A,u±;uo) 


(1 - - '^±)) (1 - - ^±)) 


Remark. Note that exp(A, ■]Vo) = e(2/A, •; vq). Hence, e and exp are equivalent up to reparametriza- 
tion. On square lattices, the discrete exponential was already considered by Ferrand m and Duffin 
[8]. The discrete exponential e on rhombic lattices was used in the work of Kenyon [l5], the first 
author, Mercat, and Suris [T], and Bucking [3j. To be comparable to their work, we use their notion 
to perform our calculations of the asymptotic behavior. In contrast, Mercat [l9] and Chelkak and 
Smirnov [5] preferred the parametrization of exp that is closer to the smooth setting. Indeed, Mercat 
remarked that the discrete exponential exp in the rhombic setting is a generalization of the formula 


exp(Ax) = 



+ 0 






to the case when the path from the origin to x consists of 0{\x\/5) straight line segments of length 5 
of any directions [El¬ 


se 











If V is fixed, exp(-,u;tio) is a rational function with poles at all edges of a shortest directed path 
connecting vq with v. It follows from Lemma IA.2I that the arguments of all poles are contained in 
an interval of length less than vr. If in addition the interior angles of parallelograms are bounded 
from below by ag, the arguments of all poles lie even in an interval of length at most tt — oq by 
Lemma fA.21 But in any case, the following motivation to call exp discrete exponential can be checked 
in a straightforward way. 

Proposition 3.1.1. For any X gG, vq G P(A), and Q G P(0); exp(A, •;uo) is discrete holomorphie 
and 

(9a exp (A,-; Wo)) {Q) = Aexp(A, Q; uq). 

3.2 Asymptotics of the discrete Green’s function 

Following the presentation of the first author, Mercat, and Suris [T], we first define a discrete loga¬ 
rithmic function on a certain branched covering A^^ of A. 

Fix Vq G F(A), and let ei,e 2 , • • •, e„ be the directed edges starting in uq, ordered according to their 
slopes. To each of these edges e we assign the real angle 9e := arg(e) G [0,27r). We assume that 
0ei < Now, define 9a+bn •= + 27r6, where a G {1,..., re} and b G'E. 

Let Ue C V (A) denote the set of all vertices to that vq can be connected by a directed path of edges 
whose arguments he in [arg(e), arg(e) + tt), where e is one of the e^. Lemma lA . 2 1 shows that the union 
of all these C/gfe) k = 1, ... ,n, covers the whole quad-graph. It follows that 

CXD 

u--= U 

m=—oo 

defines a parallelogram-graph A^^ that is a branched covering of A, branched over vq. Here, Um is 
the set of all vertices of A to that vq can be connected by a directed path of edges whose arguments 
lie in [9m, 9m + tt), equipped with the additional datum of this interval. Then, all Um+bn, h gTj and 
1 < rre < re, cover the same sector Ue^, and Um Fl Um' 7 ^ 0 if and only if |rre — m'\ < re. 

To each vertex v gV (A„g) covering a vertex u 7 ^ uq of A, we assign 9^ = arg(u — vq) mod 27r such that 
9v £ [9m, ^m + 7 r) if w G Um- Then, 9^ increases by 27r when v winds once around vq in counterclockwise 
order; and if h,!}' 7 ^ uq are adjacent vertices of A^p, [9y — 9yi[ < n. 

Note that the strip passing through an edge e separates its two endpoints from each other in the 
quad-graph A^g. For the definition of a strip, see the appendix. By construction, if we connect vq to 
some h 7 ^ uq by a shortest directed path of edges of A^^, the angles assigned to the edges lie all in 
{9y -Tr,9y + tt). 

Definition. Let vq G F(A) and let A„p be the corresponding branched covering of A. The discrete 
logarithmic function on V (A^p) is given by 

log(h;uo) ■= ^ j ^-^^^^i^U;vo)dX, 

where Cy is a collection of counterclockwise oriented loops going once around each pole of e(-,u;uo), 
V G F(A) being the projection of h G IL(A„q). The arguments of all poles shall he in {9y — 7r,9y + vr). 

It is easy to see that the real part of the discrete logarithm log(-; uq) is a well defined function on V (A). 
Divided by 27r, one actually obtains a discrete Green’s function with respect to vq. The first author. 
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Mercat, and Suris showed that this function coincides with the one of Kenyon m in the rhombic case 
[T]. Their proof can be adapted to our setting. 

Proposition 3.2.1. Let vq G K(A). The function G{-;vo) : K(A) —)• R defined by G{vo;vo) = 0 and 


a{v,v,) = if-J 
\ 


log(A) 

2A 


e(A, v; vo)dX 


is a (free) discrete Green’s function with respect to vq. Here, Gy is a collection of counterclockwise 
oriented loops going once around each pole of e{-,v,vo), where the arguments of all poles shall lie in 
(arg(u - Vo) - TT, arg(v - vq) + vr). 

Theorem 3.2.2. Assume that there are ao,qo > 0 such that a > ao and e/e' > qo for all interior 
angles a and the two side lengths e,e' of any parallelogram of A. Let vq € K(A) be fixed. 

Then, the discrete Green’s function G{-]Vo) constructed in Provosition \3.2.1\ has the following asymp¬ 
totic behavior: 


G{v,vo) 

G(v;vo) 


1 V - Vo 

— log - 

dvr J(v,vo) 


+ O (|u - Vo 


7Euler + log(2) 
27r 


+ -log|(u- 


if u and vo are of different color, 
vo)J{v, uo)| + O (|u — uo|~^) otherwise. 


Here, 7Euier denotes the Euler-Mascheroni constant. 


The proof follows the ideas of Kenyon m and Bucking [3]. Both considered just quasicrystallic 
rhombic quad-graphs. But the main difference to m is that we deform the path of integration into 
an equivalent one different from Kenyon’s, since his approach does not generalize to parallelogram- 
graphs. As Chelkak and Smirnov did for rhombic quad-graphs with bounded interior angles in [5], 
Kenyon used that two points v,v' G K(A) can be connected by a directed paths of edges such that 
the angle between each directed edge and v' — v is less than 7r/2 or the angle between the sum of 
two consecutive edges and u' — u is less than 7r/2. This is true for rhombic quad-graphs, but not for 
parallelogram-graphs. Instead, we use essentially the same deformation of the path of integration as 
Bucking did. In this paper, we will omit the straightforward calculations, the full proof can be found 
in the thesis of the second author |12j . 


Proof. The poles ei,..., e^y') of e(-, u; vo) correspond to the directed edges of a shortest path from vq 
to V. By Lemma TA. 2 1 there is a real 6o such that the angles associated to the directed edges above lie 
all in [6o, 6o + — ao]- Without loss of generality, we assume 9o = — (vr — ao)/2. By definition, 

G{v, uo) = Re J j , 

where Gy is a collection of counterclockwise oriented loops going once around each ei,... ,ek(y). The 
arguments of the poles correspond to the assigned angles. By residue theorem, we can deform Cy 
into a new path of integration G( that goes first along a circle centered at 0 with large radius R{v) 
(such that all poles lie inside this disk) in counterclockwise direction starting and ending in —R{v), 
then goes along the line segment [—i?(u), —r(u)] followed by the circle centered at 0 with small radius 
r{v) (such that all poles lie outside this disk) in clockwise direction, and finally goes the line segment 
[—R{v), —r{v)] backwards. Note that log differs by 27ri for the two integrations along [—ii(u), —r(u)]. 
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By Proposition IA.3[ there are Ei > Eq > ^ such that Ei > e > Eq for all lengths e of edges of A. In 
particular, we can choose the radii of order 

R{v) = n(|n — uol^) and r{v) = n(|n — 


By construction, k{v) < Re(n — no)/(-Eo cos(0o))) so k{v) = 0{\v — t>o|). Also 


k(v) 

i=i 


<Ma„d 


Eo 


El 


Hence, |J(n,no)| = Cl{\v — no|). 

We first look at the contributions of the circles with radii r{v) and R{v) to G{v,vo). For A —)• 0, 
(— l)^('“)e(A, V, Vo) = 1 + 0(A|t; — no|), and for A oo, e(A, v, no) = 1 + 0(A“^|n — no|). 

Thus, we get (—log(r(n))/(47 r) (l + 0(|n — no|“^)) as the contribution of the integration along 
the small circle with radius r{v). Similarly, log(i2(n))/(47r) • (l + 0{\v — no|“^)) is the contribution of 
the circle of radius R{v). 

The two integrations along [—ii(n), —r(n)] can be combined into the integral 

—r{v) 

1 r 

dvr 7 A 

-R(v) 


For |n — no| > 1 large enough, we split the integration into the three parts along 

[-R(n), -Eix^lv - no|], [-Eixy\v - no|,-], and [- y=£^=,-r{v)]. 

x/\v-vo\ VT“^o| 


We hrst consider the intermediate range A G [—Eix/\i}^^^\, —El/^^/\v^^vo\]■ Using that 


= 1 + 


|A + ep 
|A - e|2 

|e(A,n;no)| < exp 


4ARe(e) 


A2 -2ARe(e) + |e|2 

2A Re(n — vq) 


4ARe(e) 

< 1 + -r-rw 7 < exp 


(A-|e|)2 


4ARe(e) 

(A-^i)2 


< exp 


2A cos(0o)|'(^ — 
(A - i?l)2 


(A - Ei)2 

As a consequence, the contribution of the intermediate range decays faster to zero than any power of 
\v - -yol. 


Now, let A G [—Ei/y/\v — ?;ol, —r(ti)] be small. Then, 

(-l)^(’')e(A,n; uq) = exp {2XJ{v,vo) + 0{k{v)X^)) = exp {2XJ{v,vo)) (l + 0{\v - r^o|A^)) . 
Thus, the integral near the origin is equal to 

? /exp(2AJ(n,'Co)) , xx zxo l\ 2 x^J^ 

— — / I -- \-exp{2XJ{v;vo))0{\v-vo\X )jdX 

-El/y/\v-Vo\ 


fTEuler + Hd?; 


vo\ + log(2r(t>)J(t>,no))) . 
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Similarly, the integration along A E [—R{v), —Eiy^\v — uqI] gives 

TEuler - log “ ^oT^)- 

Summing up and taking the real part, we hnally get 
47rG(u;uo) = (l + (TEuier + log(2)) +log|u - uol + log |J(u;'(;o)| +0{\v - uqT^). 

□ 

Remark. Under the conditions of Theorem l3.2.2l it follows from Proposition lA.31 and Theorem l2.4.12l 
that there is exactly one discrete Green’s function having the same asymptotics as the one we con¬ 
structed above up to terms of order o{\v — 

Let us compare this result to the case of rhombi of side length one. Assume that vq E U(r). Then, 
the discrete logarithm is purely real and nonbranched on U(r) and purely imaginary and branched 
on U(r*). Thus, G{v;vo) = 0 if u E U(r*), well fitting to the fact that A splits into two discrete 
Laplacians on T and T*. Using J{v,vo) = v — vq, 

G{v; Vo) = 7 ^ (TEuier + log(2) + log \v - Vol) + 0(|v - VoT^), 

ZTT 

exactly as in the work of Bucking [3], who slightly strengthened the error term in Kenyon’s work 
m- In this paper, Kenyon showed that there is no further discrete Green’s function with asymptotics 
o{\v - Vol). 


3.3 Asymptotics of discrete Cauchy’s kernels 

Let Vo E U(A) and Qo £ C(0)- We first construct a discrete Gauchy’s kernel with respect to vo 
on U(0) that has asymptotics I2(|Q — vo|~^). Then, we construct a discrete Cauchy’s kernel Kq^ with 
respect to Qo on U(A) with asymptotics n(|v — Qo|~^)- la both cases, there are no further discrete 
Cauchy’s kernels with asymptotics o{\Q — vo|“^'^^) or o(|v — la the end of this section, we 

determine the asymptotics of OaRq^. 

Theorem 3.3.1. Let vq E U(A) and G(-;vo) be a discrete Green’s function with respect to vq. 


(i) Ky^ := 87rdAC(-;vo) is a discrete Cauchy’s kernel with respect to vq. 

(a) Assume additionally that there exist positive real numbers ao, qo such that a > ao (I'^d e/e' > go 
for all interior angles a and the two side lengths e,e' of any parallelogram of k. Suppose that 
G(-;vo) is the discrete Green’s function constructed in Provosition \3.2.1] and Ky^ the discrete 
Cauchy’s kernel given in (i). Then, 


Ky,{Q) 


1 t{Q,vo) 
Q-vo J{Q,vo) 


+ 0{\Q-vo\-^). 


(Hi) 


Under the conditions of (ii), there is exactly one discrete Cauchy’s kernel with respect to vq with 
asymptotics o{\Q — vo|“^'^^). 
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The first part is an immediate consequence of Corollary 12.4. 11 A straightforward calculation shows 
the second part, Theorem 12.4.121 the third part. 

Since we do not have discrete Green’s functions on F(<C>), we have to construct discrete Cauchy’s 
kernels on V (A) differently. To do so, we follow the original approach of Kenyon using the discrete 
exponential m that was reintroduced by Chelkak and Smirnov in [5]. Their proofs can be adapted 
to our setting. 

Proposition 3.3.2. Let Qq G K(0)- The function defined by 

0 

^Qo(^) J log(^)e(A,u;Qo)c^A = 2 J e{\,v,QQ)d\ 

Cv —{v—Qo)oo 

is a discrete Cauchy’s kernel with respect to Qq. Here, Cy is a collection of counterclockwise oriented 
loops going once around each pole of e{-, v, Qo). The arguments of all poles shall lie in {9y — vr, Oy + Ti), 
where 9y := arg(u — Qo)- 

Theorem 3.3.3. Assume that there are ao,qo > 0 such that a > oq o-nd e/e' > qq for all interior 
angles a and the two side lengths e, e' of any parallelogram of A. Let Qq G K(0) be fixed. 


(i) The discrete Cauchy’s kernel Kq^ constructed in Provosition \3.3.^ has the following asymptotics: 

1 


KqM = 




v-Qo J{v,Qo) 


(a) There is no further discrete Cauchy’s kernel with respect to Qq of asymptotics o{\v — Qo\ 

(Hi) For Kqq constructed in Proposition 1 3.3. has asymptotics 

= - (^ - jHw +° ■ 

fivj Up to two additive constants on V and T*, there is no further discrete Cauchy’s kernel with 
respect to Qo such that its discrete derivative has asymptotics o(|Q — 


The proof of the first part follows the ideas of Kenyon m- Similar to the proof of Theorem 13.2.21 
the path of integration is deformed into one different from {—{v — Qo)oo, 0] that was used by Kenyon. 
For the same reasons as before, his approach does not generalize to parallelogram-graphs. The second 
and the fourth part of the theorem are immediate consequences of Theorem 12.4.121 the third part is 
shown by a direct computation. 

Remark. Note that Kenyon and Chelkak and Smirnov proved that in the rhombic setting there is a 
unique discrete Cauchy’s kernel with asymptotics o(l) [T^[5]. 


3.4 Integer lattice 

Let us consider a locally finite planar parallelogram-graph A such that each vertex has degree four. 
This happens if and only if 0 is also a quad-graph or, equivalently, if A has the combinatorics of the 
integer lattice 


41 


















As usual, let F and F* denote the subgraphs on black and white vertices, and let us place the vertices 
of <0 at the centers of corresponding faces. Then, the vertices of 0 lie at the midpoints of diagonals. 
Since any vertex of F or F* is enclosed by a quadrilateral of F* or F, respectively, the faces of 0 are 
parallelograms by Varignon’s theorem. Thus, 0 becomes a planar parallelogram-graph as well. 

Of particular interest is the case that the two notions of discrete holomorphicity on 0, the one coming 
from 0 being the dual of A and the other coming from the quad-graph 0 itself, coincide. It is not 
hard to show that this happens only for the integer lattice of a skew coordinate system, onto which 
we restrict ourselves in the following. If ei, 62 denote two spanning vectors, 0 is a parallel shift of A 
by ei/2 -|- 62 / 2 . Furthermore, the discrete differentials on 0 seen as the dual of A coincide with the 
discrete differentials on 0 seen as a parallelogram-graph. 

Since corresponding notions coincide and 0 and A are congruent, we can skip all subscripts A and 0 in 
the definitions of discrete derivatives. Moreover, the discrete Laplacian A is now defined for functions 
on I^(A) and functions on V{()) in the same way. Due to Corollary 12.4.11 4dd = A = 4dd is now true 
on both graphs. It follows that all discrete derivatives d^f of a discrete holomorphic function / are 
discrete holomorphic themselves. Conversely, a discrete primitive exists for any discrete holomorphic 
function on a simply-connected domain by Proposition 12.2.81 

Our main interests lie in giving discrete Cauchy’s integral formulae for higher order derivatives of a 
discrete holomorphic function and determining the asymptotics of higher order discrete derivatives of 
the discrete Cauchy’s kernel constructed in Section 13.,'ll Note that due to the uniqueness statements 
in Theorems 13.3.11 and 13.3.31 both constructions yield the same discrete Cauchy’s kernel. 

Without loss of generality, we restrict our attention to functions on 14(A). For the ease of notation, 
we introduce the discrete distance D{-, ■) on 14(A) U 14(0) that is induced by the | • |oo-distance on the 
integer lattice spanned by 61 / 2 , 62 / 2 . 

Theorem 3.4.1. Let vq G 14(A), Qq := vq + 6i/2 -|- 62/2 G 14(0), let f be a discrete holomorphic 
function on 14(A) and let and Kq^ be discrete Cauchy’s kernels with respect to vq and Qq, re¬ 
spectively. Let n be a nonnegative integer and define xq := vq if n is even and xq := Qq if n is odd. 
Similarly, let x G 14(A) if n is even and x G 14(0) if n is odd. 


(i) For any counterclockwise oriented discrete contour 
x' G 14(A) U 14(0) with D{x' ,xq) < n/2, 


in the medial graph X enclosing all points 


d^fixQ) 



fd'’^Kxodz. 


(ii) If Kqq is the discrete Cauchy’s kernel constructed in Proposition 1,9. S . gl 


(-ir 


n\ 




+ 


t'{x,Qq) 


(x - ( 5 o)"'+^ (d(x, (5o)ei62)"+^ 


+ 0{\x ~ Qo\ 


—n—3\ 


where t'{x,Qq) = 1 if x and Qq or {x + ei /2 + 02/2) and Qq can be connected by a path on 14(0) 
of even length and t'{x,Qq) = —1 otherwise. 


Proof, (i) Let D be the discrete domain in X bounded by C^q. By the assumptions on Cxq, the 
discrete one-form Kx^dz vanishes on Cx^. Thus, 

f fd^Kx.dz = y" fd{d^-^Kx,) =11 d{fd{d^-^Kx,)) = 11 dfAd{d^-^Kx,) 

CxQ D D 
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by discrete Stokes’ Theorem 12.3. 11 Theorem 12.3.81 and Proposition 12.3.21 Now, / is discrete holomor- 
phic, so df A d Kxq) = dfdd^~^KxQO,(^. But since the discrete derivatives commute according 
to Corollary 12.3.31 dd^~^KxQ = d"'~^dKxQ vanishes outside Cxq, so by subsequent use of Proposi¬ 
tion [2221 

/ fd^Kx.dz = -2i{df,d^-^dKx,) = 2i{-ir{d^f,dKx,) = 2mi-ird^fixo). 

(ii) Again, the proof follows the lines of the proof of Theorem 13.2.21 In addition, we use that Propo¬ 
sition 13.1.11 shows that the discrete derivative of the discrete exponential is up to an explicit factor 
again a discrete exponential. Again, we exchange the discrete differentiation and the integration in 
the proof. □ 


A Appendix: Planar parallelogram-graphs 

The aim of this appendix is to discuss some combinatorial and geometric properties of parallelogram- 
graphs that were used in Section [3l The following notion of a strip is standard, see for example the 
book of the first author and Suris [ 2 ] . 

Definition. A strip in a planar quad-graph A is a path on its dual 0 such that two successive faces 
share an edge and the strip leaves a face in the opposite edge where it enters it. Moreover, strips are 
assumed to have maximal length, i.e., there are no strips containing it apart from itself. 

Note that a strip is uniquely determined by the edges it passes through, meaning the edges two 
successive faces share. 

Definition. For a strip S' of a parallelogram-graph A, there exists a complex vector as such that any 
(nonoriented) edge through which S passes is ± 05 . We call as a common parallel. 

as is unique up to sign; the choice of the sign induces an orientation on all edges. The parallel edges 
of the strip can be rescaled to length |a 5 | = 1, without changing the combinatorics. Hence, rhombic 
planar quad-graphs and planar parallelogram-graphs are combinatorially equivalent. Rhombic planar 
quad-graphs are characterized by the following proposition of Kenyon and Schlenker m- 

Proposition A.l. A planar quad-graph A admits a combinatorially equivalent embedding in C with 
all rhombic faces if and only if the following two conditions are satisfied: 

• No strip crosses itself or is periodic. 

• Two distinct strips cross each other at most once. 

That planar parallelogram-graphs fulfill these two conditions was already noted by Kenyon in m- 
The reason for that is that any strip S is monotone with respect to the direction ias- The coordinates 
of the endpoints of the edges parallel to as are strictly increasing or strictly decreasing if they are 
projected to ias. Whether the projections are decreasing or increasing depends on the direction in 
which the faces of S are passed through. Without loss of generality, we assume that the faces of S 
are passed through in such a way that the projections of the corresponding coordinates are strictly 
increasing. For Q £ S, let denote the semi-infinite part of S starting in the quadrilateral Q that 
passes through the faces of S in the same order. 
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As a consequence, no strip crosses itself or is periodic. Furthermore, S divides the complex plane into 
two unbounded regions, to one is as pointing and to the other —as- When a distinct strip S' crosses 
S, it enters a different region determined by S, say it goes to the one to which as is pointing. Due 
to monotonicity, the angle between ias' and as is less than 7r/2. It follows that S' cannot cross S 
another time, since it would then go to the region —as is pointing to, contradicting that the angle 
between ias' and —as is greater and not less than tt/2. 

In order to construct the discrete Green’s function and the discrete Cauchy’s kernels in Sections 13.21 
and ESI we chose a particular directed path connecting two vertices (or a face and a vertex) by 
edges of the parallelogram-graph A. This path was monotone in one direction and any angle between 
two consecutive was less than tt. The existence of such a path follows from the following lemma, 
generalizing a result of the first author, Mercat and Suris [1] to general parallelogram-graphs. The 
proof bases on the same ideas. 

Lemma A.2. Let K he a parallelogram-graph and let vq E D(A) be fixed. For a directed edge e of A 
starting in vq, consider the set Ue C D(A) of all vertices to that vq can be connected by a directed path 
of edges whose arguments lie in [arg(e), arg(e) -|- vr). 

Then, the union of all Ue, e a directed edge starting in vq, covers the whole quad-graph. 

If there is ao > 0 such that a > uq for all interior angles a of faces of A, the same statement holds 
true if [arg(e), arg(e) -|- vr) is replaced by [arg(e), arg(e) -|- vr — ao]. 

Proof. Let us rescale the edges such that all of them have length one. By this, we change neither the 
combinatorics of A nor the size of interior angles. 

For a directed edge e starting in vq, let Ufi and Ufi denote the (directed) paths on A starting in vq, 
obtained by choosing the directed edge with the least or largest argument in [arg(e), arg(e) -|- vr) (or 
[arg(e), arg(e) -|- vr — ao]) at a vertex, respectively. We hrst show that all vertices in between Ufi and 
Ufi are contained in Ue- Then, it follows that Ug is the conical sector with boundary Ufi and Ufi. 

Suppose the contrary, i.e., suppose that there is a vertex v between Ufi and Ufi to which vq cannot 
be connected by a directed path of edges whose arguments lie all in the interval [arg(e), arg(e) -|- vr) 
(or [arg(e), arg(e) -|- vr — ao]). Let the combinatorial distance between vq and v be minimal among all 
such vertices. 

In the case that interior angles of rhombi are bounded by ao from below, they are bounded from above 
by vr — ao. Hence, there is a vertex v' adjacent to v such that the argument of the directed edge v'v 
lies in [arg(e), arg(e) -|-vr — ao]. Even if interior angles are not uniformly bounded, v' can be chosen in 
such a way that the argument of v'v lies in [arg(e), arg(e) -|- vr). By construction, v' is not in C4, but 
still between Ufi and Ufi. Let us look at the strip S passing through v'v. Suppose that the common 
parallel as points from v' to v. 

If S intersects Ufi or Ufi, an edge parallel to as is contained in Ufi or Ufi, respectively. By construction, 
Vq and v' then lie on the same side of the strip S. 

If S does neither intersect Ufi nor Ufi, it is completely contained in the left half space determined by 
the oriented line vq + te, t G R, as Ufi and Ufi are. Suppose S intersects the ray vq tas, t > 0. 
Again, it follows that vq and v' he on the same side of S. 

It remains the case that S neither intersects Ufi, Ufi, nor the ray vq + tas, t > 0. Consider the 
quadrilateral area R in between the parallels vq tas, v' -\- tas and vq -\- te, v' te, t G H. By 
assumption, the semi-infinite part of S that starts with the edge v'v and then goes into R does not 
intersect an edge of R incident to vq, and by monotonicity, it does not intersect v' + tas again. Now, 


44 


A is locally finite, such that only finitely many quadrilaterals of S are inside P. Thus, S leaves P on 
the edge v' + te, t G R, and it follows that S separates vq and v. 

So in any case, S separates vq from v. Any shortest path P connecting both points has to pass through 
S. Let w be the first point of P that lies on the same side of S' as u does. Any strip passing through 
an edge on the shortest path connecting w and v on S has to intersect P as well. It follows that 
replacing all edges of P on the same side of S' as u by the path from w to v does not change its length. 
But then, v' is combinatorially nearer to vq than v, contradiction. 

Finally, we can cyclically order the directed edges starting in vq according to their slopes. Then, the 
sectors Ue are interlaced, i.e., Ue contains both U^_ and where e_,e, e+ are consecutive according 
to the cyclic order. As a consequence, the union of all these Ug covers the whole of B(A). □ 

To perform our computations in Sections 13.21 and 13.31 we needed not only that the interior angles 
were bounded, but also that the ratio of side lengths of all parallelograms were bounded. For general 
quad-graphs, we used instead boundedness of side lengths to show in Theorem 12.4.121 of Section [2] 
that biconstant functions are the only discrete harmonic function whose difference functions on 14 (F) 
and I4(r*) have asymptotics When the side lengths are bounded, the ratio of side lengths is 

trivially bounded as well. Under the assumption of bounded interior angles, the converse is true for 
parallelogram-graphs. 

Proposition A.3. Let K he a parallelogram-graph and assume that there are > 0 such that 

a > ao and e/e' > go for all interior angles a and two side lengths e, e' of any parallelogram of A. 
Then, there exist positive numbers Ei > Eq such that Ei > e > Eq for all edge lengths e. 

Proof. Let Q' G U(0) be fixed with edge lengths ei > cq and let Q G U(0) be another parallelogram 
with center x. In the following, we construct a sequence of n strips such that any two consecutive strips 
are crossing each other, the first one contains Q', the last one contains Q, and n < N := [[27r/aoJ/^l. 
Then, it follows that the side lengths of Q are bounded by Eq := and Ei := ciq^^. 

Let So be a strip containing Q'. If Q G S'o, we are done. Otherwise, we choose the common parallel 
asg in such a way that x lies in the region —as^ is pointing to. Since is monotone in the direction 
iaso and interior angles are uniformly bounded, the ray x + tasg, t > 0, intersects S in exactly one 
line segment. Let yo be the first intersection point and Qq a quadrilateral of S containing yo- 

Because A is locally finite, the line segment connecting x and yo intersects only finitely many par¬ 
allelograms. Through any such parallelogram at most two strips are passing. Thus, only a finite 
number of strips intersect this line segment. Therefore, we can choose a strip Si intersecting Sq° in a 
parallelogram Qo,i such that S'! does not contain Q and does not intersect the line segment connecting 
X and yo- Moreover, we require that Q' ^ Sq°’^. Now, choose the common parallel asi of Si in such 
a way that vr -|- arg(a 5 p) > arg(a 5 ^) > arg(a 5 p). By construction, x lies in the region — 05 ^ is pointing 
to. Note that cannot cross So a second time. 

Suppose we have already constructed the strip Sk with common parallel 05 ^, , k > 0, and x lies in the 
region — 05 ^ is pointing to. S^ shall not intersect the line segments connecting x and yo, or connecting 
X and yk-i- Moreover, assume that the semi-infinite part starting in the intersection Qk-i,k 

of Sk with Sk-i does not cross So- 

Let yk be the first intersection of the ray x + tas^, t > 0, with a quadrilateral Qk G Sk- By the same 
arguments as above, there exists a strip Sk+i intersecting n 5^'“ that does not contain Q and 

does not intersect the line segments connecting x and yo or x and yk- Choose its common parallel 
in such a way that vr -|- arg(a 5 ^) > arg(a 5 j.^J > arg(a 5 j,). By construction, x lies in the region 
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is pointing to. If the semi-infinite part starting in the intersection Qk,k+i with Sk does 

not cross Sq, we continue this procedure. 



Figure 8: Schematic picture of the proof of Proposition IA.3I 

After at most I := |'27r/aol steps, we end up with a strip Si such that intersects Sq. Indeed, 

let us suppose the contrary, that is, let us suppose that all ..., do not cross So. 

By assumption, arg(a 5 j,) + -k — ao > arg(o 5 j,^J > arg(a 5 ^) -|- ao- It follows that the first j such 
that arg(a 5 ^.) is greater or equal than arg(a 5 p) -|- 27r satishes j < |'27r/ao]. In addition, we have 
arg(a5.) < arg(a5j + 37 t - oq . 

By construction, Sj does not intersect the line segment connecting x and Vj-i- Moreover, 
cannot cross 5,_i a second time. It follows that cannot intersect the ray x -|- fa^.,, t > 0. 

Also, does neither cross So nor does it intersect the line segment connecting x and yo; so it 

does not intersect the ray x + tas^, t > 0. Thus, is contained in the cone with tip x spanned by 

and asQ (with angle less than vr). This contradicts the monotonicity of into the direction 

iasj, because the ray x + tas^, t > 0, is not contained in the interior of the cone above. 

In summary, we found a cycle of m strips Sq, Si,, Sm-i surrounding x, where m < [27r/Q!oJ + 1. 
Actually, m < [27r/Q;oJ; because the are cyclically ordered. Since only finitely many strips intersect 
the strip Sq in between Q' and Qo,ij we can assume that Q' is contained in Sq’^~^’°. 

These m strips determine a bounded region x is contained in. If Q' / we look at the semi- 

inhnite part of the strip Sq different from Sq that passes through Q' and goes into the interior of the 
bounded region above. It has to intersect one of the strips Si,..., Sm-i, say Sk- Then, Sq, ..., Sk, Sq 
or Sq, Sk,..., Sm-i, So determine a bounded region x is contained in {Q may be an element of Sq). 
Clearly, they are at most [27r/Q:oJ such strips, and Q' lies on an intersection. 

If Q ^ ^o, a strip Sq containing Q has to cross two different strips of the cycle due to local finiteness. 
In the same way as above, we can find a cycle of at most m' < [27r/aoJ strips Sq, ..., such 

that Q lies on one of the strips, say S'j^, and the intersection of Sq and is Q'. If A: < [m'/2], 

we choose the sequence of strips Sq, S'l,..., S'j^; otherwise, we take S'^,_ 2 ,..., S'j^. Any two 

consecutive strips are crossing each other, Q' is on the hrst strip, Q on the last one, and there are at 
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most |'[27r/aoJ/2l of them. 


□ 


Remark. In general, the bound [[27r/Q;oJ/2] in the proof is optimal. Indeed, consider n rays emanat¬ 
ing from 0 such that the angle between any two neighboring rays is 27r/n. In each of the n segments, 
choose the quad-graph combinatorially equivalent to the positive octant of the integer lattice that is 
spanned by two consecutive rays. For example, if n = 4, we obtain W?. Then, any strip passes through 
exactly two adjacent segments, and [n/2] is the optimal bound. 
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